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These lectures on non-baryonic dark matter matter are divided into two parts. In the 

first part, I discuss the need for non-baryonic dark matter in light of recent results in 

^ • cosmology, and I present some of the most popular candidates for non-baryonic dark matter. 

^ ■ These include neutrinos, axions, neutralinos, WIMPZILLAs, etc. In the second part, I 

^ ! overview several observational techniques that can be employed to search for WIMPs (weakly 

^ I interacting massive particles) as non-baryonic dark matter. Among these techniques, I 

\^ I discuss the direct detection of WIMP dark matter, and its indirect detection through high- 

O ' energy neutrinos, gamma-rays, positrons, etc. References cited in these lectures are intended 

2 ■ mostly for further study, and no attempt has been made to provide a comprehensive list of 

'sj- , original and recent work on the subject. 
O ■ 

1 The need for non-baryonic dark matter 

o ■ 

' We live in a time of great observational advances in cosmology, which have given us a 
ci • consistent picture of the matter and energy content of our Universe. Here matter and 
> ! energy (which special relativity tells us are equivalent) are distinguished by their different 
I dependence on the cosmic volume: matter density decreases with the inverse of the volume, 
' while energy density remains (approximately) constant. 
■ — ' Nothing is known about the nature of the energy component, which goes under the 

name of dark energy. Of the matter component, less than 2% is luminous, and no more 
than 20% is made of ordinary matter like protons, neutrons, and electrons. The rest of 
the matter component, more than 80% of the matter, is of an unknown form which we call 
non-baryonic. Finding the nature of non-baryonic matter is referred to as the non-baryonic 
dark matter problem. 

A summary of the current measurements of the matter density and the energy 
density Q,\ are shown in Figure 1 (adapted from Verde et al. . 2002f l. Both are in units of 



the critical density pcrit = S-ffg/ (SttG), where G is the Newton's gravitational constant, and 
Hq is the present value of the Hubble constant. Three types of observations - supernova 
measurements of the recent expansion history of the Universe, cosmic microwave background 
measurements of the degree of spatial flatness, and measurements of the amount of matter 
in galaxy structures obtained through big galaxy redshift surveys - agree with each other 
in a region around the best current values of the matter and energy densities Q^a — 0.27 
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Figure 1: The concordance cosmology and the need for non-baryonic dark matter. Current 
cosmological measurements of the matter density and energy density Qa give the value 
marked with a cross at Q^i — 0.27, Q\ ~ 0.73. The baryon density does not exceed 
0.05 (black vertica l band). The rest of the matter is non-baryonic. (Figure adapted from 



Verde et ahLliooi^ 



and — 0.73 (cross in Figure 1). Measurements of the baryon density in the Universe 
using the cosmic microwave background spectrum and primordial nucleosynthesis constrain 
the baryon density fib to a value less than ~ 0.05 (black vertical band in Figure 1). The 
difference flm — fib — 0.22 must be in the form of non-baryonic dark matter.^ 

A precise determination of the cosmological density parameters is able to give the matter 
and en ergy densities i i i phys ical units. For example, in units of 1.879 x 10~^^ g/cm^ = 18.79 
yg/m^, Spergel et al. ( 2003[ ) have determined a total matter density 
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^The red vertical band labeled 'stars' in Figure 1 shows the density of lu minous matter, corrected for 
all expected dim stars and gas (see, e.g.. iFukugita. Hogan. fc PeeblesL Il998j) . The difference between the 
amount of luminous matter and the amount of baryons constitutes the dark baryon problem, which will not 
be addressed here (see, e.g., 



fi^/i^ < 0.0076 

is in the form of neutrinos (to 95% confidence level), 

fib/i^ = 0.0224 ± 0.0009 
is in the form of baryons (protons and nucleons in cosmological parlance), and 
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is in the form of cold dark matter (CDM), a non-baryonic component whose nature we are 
still trying to uncover. Some ideas of what it may be are presented in the next Section. 

2 Popular candidates for non-baryonic dark matter 

A new kind of elementary particle has been the dominant (exclusive?) candidate for non- 
baryonic dark matter. 

A major classification of non-baryonic dark matter is based on its temperature at the time 
of galaxy formation, which occurs at a photon temperature of about 1 keV. Hot dark matter 
was relativistic at the time of galaxy formation, and as a consequence hindered the formation 
of the smallest objects by streaming out of the forming structures. An example of a hot 
dark matter particle is a light neutrino, much lighter than ~keV. Cold dark matter was non- 
relativistic when galaxies formed, and thus was able to collapse effectively under the action 
of gravity because of its negligible pressure. Examples of cold dark matter particles are 
neutralinos, axions, WIMPZILLAs, solitons (B-balls and Q-balls), etc. Warm dark matter 
was semi-relativistic at the time of galaxy formation, and is therefore an intermediate case 
between hot and cold dark matter. Two examples of warm dark matter are keV-mass sterile 
neutrinos and gravitinos. 

Another important classification of particle dark matter rests upon its production mech- 
anism. Particles that were in thermal equilibrium in the early Universe, like neutrinos, 
neutralinos, and most other WIMPs (weakly interacting massive particles), are called ther- 
mal relics. Particles which were produced by a non-thermal mechanism and that never had 
the chance of reaching thermal equilibrium in the early Universe are called non-thermal 
relics. There are several examples of non-thermal relics: axions emitted by cosmic strings, 
solitons produced in phase transitions, WIMPZILLAs produced gravitationally at the end 
of inflation, etc. 

For the sake of presentation, we find still another classification useful. We will divide 
candidates for particle dark matter into three categories: Type la. Type lb, and Type II 
(following a common practice in superconductors and supernovas). Type la candidates are 
those known to exist, foremost among them are the neutrinos. Type lb candidates are 
candidates which are still undiscovered but are 'well-motivated.' By this we mean that (1) 
they have been proposed to solve genuine particle physics problems, a priori unrelated to 
dark matter, and (2) they have interactions and masses specified within a well-defined (and 
consistent) particle physics model. We are aware of the arbitrariness of this classification, 
and reserve the honor of belonging to the Type lb category only to a sterile neutrino, the 
axion, and the lightest supersymmetric particle (which may be a neutralino, a gravitino, or a 
sneutrino) Finally, Type II candidates are all other candidates, some of which are examples 
of maybe fruitful ideas, such as WIMPZILLAs, solitons (B-balls, Q-balls), dark matter from 
extra-dimensions, self-interacting dark matter, string-inspired dark matter, string-perspired 
dark matter, etc. It goes without saying that a candidate may move up from Type II to 
Type lb and even to Type la as our understanding of particle physics models progresses. 

We now examine some of the current candidates. 



3 



2.1 Type la: candidates that exist 

Dark matter candidates that are known to exist in Nature have an obvious advantage 
over candidates that have not been detected. The chief particles in this category are the 
neutrinos. 

There are three known 'flavors' of neutrinos: the electron neutrino Vf.-, fhe muon neutrino 
z/^, and the tau neutrino Ur- They are so named because they are produced or destroyed in 
concomitance with the electron, the muon, and the tau lepton, respectively. 

"If neutrinos had a mass, they would be a good candidate for the dark matter," said 
Steven Hawking. We now know that neutrinos, or at least some of the neutrinos, do have a 
mass. This was discovered indirectly through the observation of neutrino flavor oscillations, 
i.e. the spontaneous conversion of one neutrino flavor into another as a neutrino propagates 
from point to point. The connection between flavor oscillations and neutrino masses can be 
seen as follows. 

Consider for simplicity two flavors of neutrinos instead of three, and z/^, say. Weak 
interactions produce the flavor eigenstates |z/e) and which are associated with their 
respective charged leptons. However, these flavor eigenstates are not energy eigenstates. 
Let and \v2) denote the two energy eigenstates for the two-flavor system, with energies 
El and E2 respectively. Then the flavor and the energy eigenstates are connected by a 
unitary transformation, 

Ve) = COS 6 lui) + sin^ W2), 

(5) 

Ui^i) = — sin 9 \ iyi) + cos 6* lz/2). 
Imagine that at time t = we produce a i^e, so that the initial wave function is 

= We). (6) 

After a time t, the wave function evolves according to the system Hamiltonian H as (we use 
natural units h = c = 1) 

= e-'^'\m) = e-''^Ve). (7) 

To see the evolution explicitly, we expand lu^) into energy eigenstates, and obtain 

\ip{t)) = cos^e-^^^Vi) + sin^e-^^^Va). (8) 

We can now ask what is the probability of observing the neutrino in the state |z/^) after 
a time t, i.e. of observing a neutrino with flavor z/^ instead of the initial Ug. According to 
standard rules of quantum mechanics, this probability is 

FToh{u,^u^) = \{u^mt))\^ (9) 

= |cos^e-*^i*(z/^|z/i) + sin^e-*^2*(z/^|z/2)f (10) 

= I -cosesin^e"*^i* + sin^cos^e"^-^2*f (11) 

= sm^29 sin^fK^a - Ei)t] , (12) 
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where we have used Eq. (0) and = (^'2|^'2) = 1, (^^11^2) = 0. For free relativistic 

neutrinos, with momentum p much larger than their mass m, we have 



E = _|_ ^2 ~ p -|- 



m 
2p 



and 



Hence 



E2 — El 



rrio 



2p 



Prob(i/e 



sin 29 sin 



ml)t 



4p 



(13) 



(14) 



(15) 



This equation shows that the probabihty of conversion from flavor z/g to flavor z/^ oscillates 
in time with a frequency proportional to the difference of the squares of the neutrino masses 
Am^2 = ''^2 '^i- The observation of neutrino flavor oscillations therefore implies that 
neutrino masses differ from each other, and in particular that at least one of them is different 
from zero. 

Neutrino oscillations have up to now been detected in two systems. Atmospheric muon 
neutrinos, which originate from the collision o f cosmic rays \yith th e Earth atmosphere, have 
been observed to oscillate into tau neutrinos ( Fukuda et al. . 1998| ). 



^23 



3 X 10" 



eVl 



(16) 



Solar neutri nos, produced in th e nuclear reactions that make the Sun shine, also show 
oscillations ( Ahmad et al. . 20021 ). 



or Ur, 



Am?2 ~ 7 X 10-^ eVl 



:i7) 



These results can be used to set a lower limit on the mass of the heaviest neutrino. 
Indeed, the mass of the heaviest neutrino must be greater than or equal to the square root 
of the largest mass-squared difference (just take the mass of the other neutrino to vanish). 
This gives the lower limit 

mass of heaviest neutrino ^ 0.05 eV. (18) 
Upper limits on neutrino masses come from laborato ry experiments, such as tritium de- 



cay a nd high-energy accelerator experiments, and are fsee lReview of Particle Phvsics. Hagiwara et al 
2002[) 

mi < 2.8 eV, < 190 keV, mg < 18.2 MeV. (19) 

However, the small mass differences implied by Eqs. (|16|) and (jl7p imply that the smallest 
of the three upper limits applies to all three active neutrino masses. Thus we have 



m,- < 2.8 eV 



1,2,3). 



(20) 



It follows from this mass constraint that reactions such as z/gi^e ^ e+e in the hot early 
universe were able to keep standard- mo del neutrinos in thermal equilibrium. The neutrino 
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density then follows from a computation of the neutrino number density (see Section 2.2). 
The result is 

(21) 

^90eV' ^ ' 

1=1 

where Qi = 1 for a neutrino which is its own antiparticle (Majorana neutrino) and Qi = 2 
for a neutrino which is not its own antiparticle (Dirac neutrino). 

We already mentioned in Section 1 that cosmology provides an upper limit on the neu- 
trino density Q^h'^. This translates into a cosmological upper limit on the neutrino mass 
using Eq. (j2H) . The cosmological limit is strictly speaking on the mass density in relativistic 
particles at the time of galaxy formation. An excessive amount of relativistic particles when 
galaxies form, i.e. of particles with mass m <^ keV, would erase too much structure at the 
smallest scales. A combination of cosmic microwave background measurements, galaxy clus- 
tering measurements, and o bservations of the Lyman-a forest gives the upper limit quoted 
before (|SDergel et al.l . 120031 ) 



n^h^ < 0.0076 (95% C.L.). (22) 

Eq. (j2H) then gives 

giuii + 5(27712 + gsm^ < 0.7 eV. (23) 

On the other hand, Eq. ()2H) can be used in conjunction with inequality ()18|) to obtain a 
lower bound on the cosmological density in neutrinos. Taking only one massive Majorana 
flavor, 

n^h"^ ^ 0.0006. (24) 

Thus neutrinos are definitely a form of dark matter, although perhaps a minor component 
of it. 

The results for the known neutrinos as dark matter can be summarized by the constraints 

0.05 eV < mi + m2 + m3 < 0.7 eV, (25) 
0.0006 < n^h^ < 0.0076, (26) 

where the constraint in Eq. (|23|) has somewhat been relaxed by taking = 1. 

The upper limit on Q^h'^ forbids currently known neutrinos from being the major con- 
stituents of dark matter. Moreover, since they are light and relativistic at the time of galaxy 
formation, the three neutrinos known to exist are hot, not cold, dark matter. 

The three active neutrinos are our only known particle candidates for non-baryonic dark 
matter. Since they fail to be cold dark matter, we are lead to consider hypothetical particles. 



2.2 Type lb: 'well- motivated' candidates 

We will discuss three cold dark matter candidates which are 'well- motivated', i.e. that have 
been proposed to solve problems in principle unrelated to dark matter and whose proper- 
ties can be computed within a well-defined particle physics model. The three candidates 
we discuss are: (1) a heavy active neutrino with standard model interactions, (2) the neu- 
tralino in the minimal supersymmetric standard model, and (3) the axion. Examples of 
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other candidates that can be in cluded in this category are a sterile neutrino (see, e.g., 



lAbazaiiaii. Fuller, fc Jr^atel l^uuil l and otner supersymm etric particles sucn as tne gravitip 
see, e g.. lEUis et al. . 1984[ ) and the sneutrino (see, e.g.. Hall. Moroi. fc MuravamaL Il998[ ). 



ate|, l200lh and other supersymm etric particles such as the gravitin o 



The first two candidates we discuss belong to a general class called weakly interacting 
massive particles (WIMPs).^ WIMPs that were in thermal equilibrium in the early universe 
(thermal WIMPs) are particularly interesting. Their cosmological density is naturally of the 
right order of magnitude when their interaction cross section is of the order of a weak cross 
section. This also makes them detectable in the laboratory, as we will see later. In the early 
Universe, annihilation reactions that convert WIMPs into standard model particles were 
initially in equilibrium with their opposite reactions. As the universe expanded, and the 
temperature became smaller than the WIMP mass, the gas of WIMPs, still in equilibrium, 
diluted faster than the gas of standard model particles. This occurred because the equilib- 
rium number density of non-relativistic particles is suppressed by a Boltzmann factor e~™/-^ 
with respect to the number density of relativistic particles. After a while, WIMPs became 
so rare that the WIMP annihilation reactions could no longer occur (chemical decoupling), 
and from then on the number density of WIMPs decreased inversely with volume (or in 
other words, the number of WIMPs per comoving volume remained constant). Chemical 
decoupling occurs approximately when the WIMP annihilation rate Fann = {o'annv)n be- 
came smaller than the universe expansion rate H. Here (Tann is the WIMP annihilation 
cross section, v is the relative velocity of the annihilating WIMPs, n is the WIMP number 
density, and the angle brackets denote an average over the WIMP thermal distribution. 
Using Friedmann's equation to find the expansion rate H gives 

ah^ « ' (27) 

for the relic density of a thermal WIMP. An important property of this equation is that 
smaller annihilation cross sections correspond to larger relic densities ("The weakest wins.") 
This can be understood from the fact that WIMPs with stronger interactions remain in 
chemical equilibrium for a longer time, and hence decouple when the universe is colder, 
wherefore their density is further suppressed by a smaller Boltzmann factor. Figure 2 
illustrates this relationship. 

It must be remarked here that in the non-relativistic limit w — ^ 0, the product o"ann^ 
tends to a constant, because the annihilation cross section a^nn diverges as 1/u as w 0. 
This is analogous to what happens for the scattering cross section of thermal neutrons. 



Heavy neutrino 

The WIMP par excellence is a heavy neutrino. The example we consider is a thermal 
Dirac neutrino u of the fourth generation with Standard Model interactions and no lepton 
asymmetry. Figure 3 summarizes its relic density as a function of mass. Also shown in the 
Figure are the current constraints from accelerator experiments and dark matter searches. 

A neutrino lighter than ~ 1 MeV decouples while relativistic. If it is so light to be 
still relativistic today (m^, ^0.1 meV), its relic density is p^, = 77r^T^/120. If it became 

^Notice that according to the Merriam- Webster Dictionary of the Enghsh Language, a wimp is a weak, 
cowardly, or ineffectual person. 
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non-relativistic after decoupling, its relic density is determined by its equilibrium number 
density as = m^3C{3)Tf/27T^. Here = {3/11)^^'%, where = 2.725 ± 0.002K is the 
cosmic microwave background temperature. (We use natural units, c = h = 1.) 

A neutrino heavier than ~ 1 MeV decouples while non-relativistic. Its relic density is 
determined by its annihilation cross section, as for a general WIMP (see Eq. (12 7p ). The 
shape of the relic density curve in Figure 3 is a reflection of the behavior of the annihilation 
cross section. The latter is dominated by the Z-boson resonance at m,^ ~ mz/'i. This 
resonant annihilation gives the characteristic V shape to the relic density curve. Above 
my ~ 100 GeV, new annihilation channels open up, namely the annihilation of two neutrinos 
into two Z- or W-bosons. The new channels increase the annihilation cross section and thus 
lower the neutrino relic density. Soon, however, the perturbative expansion of the cross 
section in powers of the (Yukawa) coupling constant becomes untrustworthy (the question 
mark in Figure 3). An alternative unitarity argument limits the Dirac neutrino relic density 
to the dashed curve on the right in the Figure. Neutrinos heavier than 10 TeV 'overdose' 
the universe, i.e. have a relic density that corresponds to a universe which is too young. 

The 'dark matter' band in Figure 3 indicates where the neutrino is a good dark matter 
candidate (the band is actually quite generous in light of the most recent measurements of 
Qh?). A thermal Dirac neutrino is a good dark matter candidate when its mass is around 
few eV, a few GeV or possibly a TeV. For masses smaller than about an eV and between 
~10 GeV and ~100 GeV, it is an underabundant relic from the Big Bang, too dilute to 
be a major component of the dark matter but nevertheless a cosmological relic. For other 
masses, it is cosmologically excluded. 

Dark matter neutrinos with a mass around 1 eV would be relativistic at the time of 
galaxy formation (~ keV), and would thus be part of hot dark matter. From the bounds 
on hot dark matter in the preceding Section, however, they cannot be a major component 
of the dark matter in the Universe. 

Figure 2: Evolution of a typical 
WIMP number density in the early 
universe. The number of WIMPs in 
a volume expanding with the uni- 
verse (comoving density) first de- 
creases exponentially due the Boltz- 
mann factor e~'^l'^ and then 'freezes 
out' to a constant value when the 
WIMP annihilation reactions cannot 
maintain chemical equilibrium be- 
tween WIMPs and standard model 
particles. In the figure, {av) is 
the thermally averaged annihilation 
cross section times relative veloc- 
ity. WIMPs with larger annihilation 
cross section end up with smaller 
densities. 
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Figure 3: Neutrinos as 
dark matter. Relic density 
of a thermal Dirac neu- 
trino with standard-model 
interactions, together with 
current constraints from 
cosmology, accelerators 
(LEP), and dark matter 
searches. See text for 
explanations. (The 'dark 
matter' band is quite 
generous in light of the 
WMAP measurements.) 
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Neutrinos can be cold dark matter if their masses are around few GeV or a TeV. However, 
fourth-generation heavy neutrinos lighter than 45 GeV are excluded by the measurement of 
the Z-boson decay width at the Large Electron-Positron collider at CERN. Moreover, direct 
searches for WIMP dark matter in our galaxy exclude Dirac neutrinos heavier than ~0.5 
GeV as the dominant component of the galactic dark halo (see Figure 3). Thus although 
heavy Dirac neutrinos could still be a tiny part of the halo dark matter, they cannot solve 
the cold dark matter problem. 

We need another non-baryonic candidate for cold dark matter. 

Neutralino 

The WIMP par default is the lightest neutralino x?? or sometimes simply x, which is often 
the lightest super symmetric particle in supersymmetric extensions of the Standard Model 
of particle physics. Supersymmetry is a new symmetry of space-time that has been discov- 
ered in the process of unifying the fundamental forces of nature (electroweak, strong, and 
gravitational). Supersymmetry also helps in stabilizing the masses of fundamental scalar 
particles in the theory, such as the Higgs boson, a problem, called the hierarchy problem, 
which basically consists in explaining why gravity is so much weaker than the other forces. 

Of importance for cosmology is the fact that supersymmetry requires the existence of a 
new particle for each particle in the Standard Model. These supersymmetric partners differ 
by half a unit of spin, and come under the names of sleptons (partners of the leptons), squarks 
(partners of the quarks), gauginos (partners of the gauge bosons) and higgsinos (partners 
of the Higgs bosons). Sleptons and squarks have spin 0, and gauginos and higgsinos have 
spin |. 

If supersymmetry would be an explicit symmetry of nature, superpartners would have the 
same mass as their corresponding Standard Model particle. However, no Standard Model 
particle has a superpartner of the same mass. It is therefore assumed that supersymmetry, 
much as the weak symmetry, is broken. Superpartners can then be much heavier than their 
normal counterparts, explaining why they have not been detected so far. However, the 
mechanism of supersymmetry breaking is not completely understood, and in practice it is 
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implemented in the model by a set of supersymmetry-breaking parameters that govern the 
values of the superpartners masses (the superpartners couplings are fixed by supersymme- 
try). 

The scenario with the minimum number of particles is called the minimal supersymmetric 
standard model or MSSM. The MSSM has 106 parameters beyond those in the Standard 
Model: 102 supersymmetry-breaking parameters, 1 complex supersymmetric parameter 
and 1 c omplex electroweak symmetry-breaking paramete r tan/3 (see, e.g., the article by 
Haber in lReview of Particle Phvsics. Hagiwara et al.l . l2002[ ). Since it is cumbersome to work 
with so many parameters, in practice phenomenological studies consider simplified scenarios 
with a drastically reduced number of parameters. The most studied case (not necessarily the 
one Nature has chosen) is minimal supergravity, which reduces the number of parameters 
to five: three real mass parameters at the Grand Unification scale (the scalar mass mo, the 
scalar trilinear coupling Aq, and the gaugino mass mi/2) and two real parameters at the 
weak scale (the ratio of Higgs expectation values tan/5 and the sign of the fi parameter). 
Other scenarios are possible and are considered in the literature. Of relevance to dark matter 
studies is, for example, a class of models with seven parameters specified at the weak scale: 
II, tan/?, the gaugino mass parameter M2, the mass of the CP-odd Higgs boson, the 
sfermion mass parameter m, the bottom and top quark t riline ar couplings Ah an d At. See 
the reviews by Jungman. Kamionkowski &: Griest ( 1996| ) and Bergstrom ( 2000l ) for more 
details. 

It was realized long ago by Goldberg ( 19831 ) and Ellis et all ()l984l ) that the lightest 
superposition of the neutral gauginos and the neutral higgsinos (which having the same 
quantum numbers mix together) is an excellent dark matter candidate. It is often the 
lightest supersymmetric particle, it is stable under the requirement that superpartners are 
only produced or destroyed in pairs (called R-parity conservation), it is weakly interacting, 
as dictated by supersymmetry, and it is massive. It is therefore a genuine WIMP, and it is 
among the most studied of the dark matter candidates. Its name is the lightest neutralino. 

Several calculations exist of the density of the lightest neutr a lino. An example is given 
in Figure 4, which reproduces a figure from Edsio fc Gondolol (1997), updated with the 
WMAP value of the cold dark matter density. This figure was obtained in a scenario with 
seven supersymmetric parameters at the weak scale. The relic density is not fixed once 
the neutralino mass is given, as was the case for a Dirac neutrino, because the neutralino 
annihilation cross section depends on the masses and composition of many other supersym- 
metric particles, thus ultimately on all supersymmetric parameters. Therefore the density 
in Figure 4 is not a single-valued function of the neutralino mass, and the plot must be 
obtained through an extended computer scan in the seven-dimensional parameter space. ^ 
It is clear from Figure 4 that it is possible to choose the values of the supersymmetric 
parameters in a way that the neutralino relic density satisfies the current determination 
of nh^. Although it may seem ridiculous to claim that the neutralino is naturally a good 
dark matter candidate, let us notice that the neutralino relic density in Figure 4 and the 
neutrino relic density in Figure 3 have a similar range of variation. It is the precision of the 
cosmological measurements that make us think otherwise. 

■^We may ask if there can be points in the empty regions, or in more general terms, what is the meaning 
of the density of points in F igure 4, and in similar figm^es in Section 3. For a discussion on this, see 
iBergstrom fc Gondolol l)l996t) . 
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Figure 4: Relic density of the 
lightest neutralino func- 
tion of its mass. For each 
mass, several density values 
are possible depending on the 
other supersymmetric param- 
eters (seven in total in the 
scenario plotted). The color 
code shows the neutralino com- 
position (gaugino, higgsino or 
mixed). The gray horizon- 
tal line is the current er- 
ror band in the WMAP mea- 
surement of the cosmologi- 
cal cold dark matter den- 
sitv- (Figure adapte d from 
Edsiofc Gondokl FmOTll 



Reversing the argument, the precision of the cosmological measurements can be used 
to select the regions of supersymmetric parameter space where the lightest neutralino is 
cold dark matter. With the current precision of cosmological measurements, these are 
very thin regions in supersymmetric parameter space. For this approach to be carried 
out properly, the theoretical calculation of the neutralino relic density should match the 
precision of the cosmological data. The latter is currently 7%, as can be gathered from 
Eq. (jH), and is expected to improve to about 1% before the end of the decade with the 
launch of the Planck mission. A calculation of the neutralino relic density good to 1% 
now exists, and is available in a computer package called DarkSUSY ( Gondolo et al. . 200(1 



120021 ). The 1% precision refers to the calculation of the relic density starting from the 
supersymmetric parameters at the weak scale. The connection with the parameters at the 
Grand Unification scale, which is vital for minimal supergravity, introduces instead large 
errors in im portant regions of parameter spac e, errors that some authors estimate to be as 



errors m im portant regions oi parameter spac 
big as 50% flAllanach. Kraml. fc Por^ \200± 



Given the importance of this calculation, we give now a rapid survey of the ingredients 
needed to achieve a precision of 1%. Firstly, the equation governing the evolution of the 
neutralino number density n, 

h + 3Hn = - (ffannt^) - <) , (28) 

has to be solved numerically. The equation being 'stiff' (i.e. its difference equation being sta- 
ble only for unreasonably small stepsizes), a special numerical method should be used. The 
thermal average of the annihilation cross section (as^nnv) at temperature T should be com- 
puted relativistically, since the typical speed of neutralinos at decoupling is of the order of the 



speed of light, f ~ c/3. For this purpose, we can use the expression in lGondolo fc Gelmini 
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J. Edsjo, M. Schelke. P. Ullio and P. Gondola, 2003 



Figure 5: Illustration 
of the power of WMAP 
constraints on the 
minimal supergravity 
parameter space. The 
figure shows a slice 
in mo and mi/2 with 
/i < at tan/9 = 30 and 
A = 0. The WMAP 
constraint is a very thin 
(grey) line that approxi- 
mately follows the edges 
of the allowed region. 

(Figure adapted from 

Edsi5. Schelke. Ullio. fc Gondol^l ?^ 
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( 199ll) : 

^ ' m^T[K2{m/T)f ^ ' 

where W{p) is the annihilation rate per unit volume and unit time (a relativistic invariant), 
s = 4(m^ + P^) is the center-of-mass energy squared, and Ki, K2 are modified Bessel 
functions. In passing, let us remark that the common method of expanding in powers of v"^, 
CannV = « + bv'^ + " " " , and then taking the thermal average to give (o"ann'?^) = a + H — ■ is 
unrehable, since it gives rise to negative (o"anni'), and thus negative Qh"^, near resonances and 
thresholds. These are nowadays the most important regions of parameter space. Finally, 
an essential ingredient in the calculation of the neutralino relic density is the inclusion 
of coannihilation processes. These are processes that deplete the number of neutralinos 
through a chain of reactions, and occur when another supersymmetric particle is close in 
mass to the lightest neutralino (Am ~ T). In this case, scattering of the neutralino off 
a particle in the thermal 'soup' can convert the neutralino into the other supersymmetric 
particle close in mass, given that the energy barrier that would otherwise have prevented it 
(i.e. the mass difference) is easily overcome. The supersymmetric particle participating in 
the coannihilation may then decay and/or react with other particles and eventually effect 
the disappearance of neutralinos. We give two examples. Coannihilation with charginos 
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(partners of the charged gauge and Higgs bosons) may proceed via, for instance, 



X2 



X2 xldu 



xlxl 



(30) 



(subscripts on superpartner names indicate particles with different masses). Coannihilation 
with tau sleptons f may instead involve the processes 



X\r 



rXi 



tW- 



(31) 



Coannihilations were first includ ed in the study of near-degenerate hea vy neutrinos by 
Binetruv. Girardi. &: Salati (1984) and were brought to general attention by Griest fc Seckell 
()l991 ). The current state-of-the-art treatment of neutralino c oannihilations, wh i ch inv olves 
several thousands of processes, i s conta ined in the work by Edsio fc Gondolol (1997) and 
Edsio. Schelke. UUio. fc Gondolol \2m± . 

To illustrate the power of the cosmological precision measurements in selecting regions 
of supers ymmetric parameter space. Figure 5 sh ows the WMAP constraint in one of the 
figures in lEdsio. Schelke. UUio. fc Gondolol (I2nn3t) . The constraint is as a very thin line that 
approximately follows the edges of the allowed region in this slice of parameter space (the 
'^o^"^i/2 plane with /U < 0, tan/3 = 30 and = 0). 

Cosmological constraints on supersymmetric models are very powerful, and may even 
serve as a guidance in searching for supersymmetry. This partially justifies the extensive 
literature on the subject. The neutralino as dark matter is certainly 'fashionable.' 



Axion 

Our third and last example of a 'well-motivated ' cold dark matter candidate is the axion. 

Axions were suggested by IPeccei fc QuinnI (jl977l ) to solve the so-called "strong CP 
problem". Out of the vacuum structure of Quantum Chromodynamics there arises a 
large CP-violating phase, which is at variance with stringent measurements of the elec- 
tric dipole moment of the neutron, for example. A possible solution to this problem 
is that the CP-violating phase is the vacuum expectation value of a new field, the ax- 
ion, which relaxes dynamically to a very small value. The original axion model of Pec- 
cei and Quinn is today experimentally ruled out, but other axion mod els based on the 
same idea have been proposed. Among them are the invisible axions of Kim ' ("1979 ^ and 
Shifnaan. Va inshtein. fc Zakh arovl (jl980l ) (KSVZ axion) and of Dine. Fischler. fc Srednicki 
()l98ll ) and IZhitnitskvi (fl98Q ) (DFSZ). They differ in the strength of the axion couplings to 
matter and radiation. 

In a cosmological context, axions, contrary to neutrinos and neutralinos, are generally 
produced non-thermally (although thermal axion production is sometimes considered, as 
are non-thermal neutrino and neutralino productions). The two main mechanisms for non- 
thermal axion production are vacuum alignment and emission from cosmic strings. In the 
vacuum alignment mechanism, a potential is generated for the axion field at the chiral 
symmetry breaking, and the axion field, which can in principle be at any point in this 
potential, starts moving toward the minimum of the potential and then oscillates around it. 
Quantum-mechanically, the field oscillations correspond to the generation of axion particles. 
In the other main non-thermal mechanism for axion production, axions are emitted in 
the wiggling or decay of cosmic strings. In both cases, axions are produced with small 
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Figure 6: Laboratory, astrophysi- 
cal, and cosmological constraints on 
the axion mass m^. The inflation 
scenario and the string scenario are 
referred in the text as the vacuum 
alignment scenario and the string 
emission scenario, respectively. 
/a is the axion decay constant, 
which is inversely related to rriA- 
The axion is a good dark matter 
candidate for 1 /leV ^ <^ 
1 meV. (Figure from Raffelt in 
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momentum, ^ keV, and thus they are cold dark matter despite having tiny masses, between 
1 neV and 1 meV. This is in fact the range of masses in which axions are good dark matter 
candidates. Figure 6 shows the current constraints on the axion mass from laboratory, 
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Figure 7: Experimental constraints 
on the density of axions in the 
galactic halo near the Sun 
function of the axion mass (upper 
scale) and cavity frequency (lower 
scale) . The regions above the curves 
marked 'DFSZ' and 'KSVZ' are ex- 
cluded fro the respective axion mod- 
els. The currently accepted value 
for the local dark halo density is 
0.45 GeV/cm^, which is approx- 
imately the extension of the ex- 
cluded regio n for the KSVZ axion . 
(Figure from Asztalos et al. . 2004 ) 
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astrophysical, and cosmological data. 

Searches for axions as galactic dark matter rely on the coupling of axions to two photons. 
An incoming galactic axion can become a photon in the magnetic field in a resonant cavity. 
For this to happen, the characteristic frequency of the cavity has to match the axion mass. 
Since the latter is unknown, searches for galactic axions use tunable cavities, and scan 
over the cavity frequency, a time-consuming process. The U.S. axion search at Livermore is 
currently exploring a wide range of interesting axion masses, a nd has put some cons traint on 
the KSVZ axion as a dominant component of the galactic halo (jAsztalos et al.Ll2004 ). Figure 
7 shows the constraints on the local galactic density in axions as a function of the axion 
mass. KSVZ axions with mass in the range 1.91-3.34 /xeV cannot be the main component 
of galactic dark matter. The Livermore search is still continuing to a larger range of axion 
masses. It is fair to say that axion dark matter is either about to be detected or about to 
be ruled out. 



2.3 Type II: other candidates 



In the Type II category we put all hypothetical cold dark matter candidates that are neither 
Type la nor Type lb. Some of these candidates have been proposed for no other reason 
than to solve the dark matter problem. Others are examples of beautiful ideas and clever 
mechanisms that can provide good possibilities for non-baryonic dark matter, but in some 
way or another lack the completeness of the theoretical particle physics models of Types la 
and lb. Although Type II candidates are not studied as deeply as others, it may well be 
that eventually the question of the nature of cold dark matter might find its answer among 
them. 

Below we present the idea of self-interacting dark matter and gravitationally-produced 
WIMPZILLAs. Other interesting candidates have be en proposed recently in models with 
extra dir nensions, such as Kaluza-Klein dark rn atter ( Cheng. Feng, fc Matchev . 2002[ ) and 
branons ( Cembranos. Dobado. fc Marotol 20031 ). 

There are several ways in which one may be able to come up with an ad hoc candidate 
for non-baryonic cold dark matter. A humorous flowchart on how to do this was put together 
aroun d 1986 by a group of graduate students at Princeton (jLauer. Statler. Rvden. fc Weinberg . 
1986t ). The flowchart involves multiple options, and one possibility runs as follows. "A new 
particle is envisioned which is cooked up just to make everything OK but violates federal law 
. . . still, it doesn't prevent a paper being writ ten by Spergel ..." We have now the proper 
setting to introduce a candidate suggested bv lSpergel fc Steinhardt self-interacting 
dark matter. 



Self-interacting dark matter 

The idea behind the introduction of self-interacting dark matter is to find a solution to the 
cusp and satellite problems of standard cold dark matter scenarios. These two problems are 
in effect discrepancies between observations and results of simulations of structure formation 
on the galactic scale. Namely, in the cusp problem, numerical simulations predict a dark 
matter density profile which increases toward the center of a galaxy like a power law r~'^ 
with 7 ~ 1 or higher. This sharp density increase is called a cusp. On the other hand, 
kinematical and dynamical determinations of the dark matter profile in the central regions 
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of galaxies, especially of low surface brightness galaxies, tend not to show such a sharp 
increase but rather a constant density core. The observational situation is still rather 
confused, with some galaxies profiles being compatible with a cusp and others with a core. 
The theoretical situation is also not fully delineated, with higher resolution simulations 
showing a dependence of the slope 7 on the mass of the galaxy. Although many ideas have 
been proposed for the resolution of the cusp problem, it is still not completely resolved. 

Perhaps connected with the cusp problem is the satellite problem, which is a mismatch 
between the observed and the simulated numbers of satellites in a galaxy halo. Too many 
satellites are predicted by the simulations. In reality, observations can detect only the 
luminous satellites while simulations contain all satellites, including the dark ones. It may 
be that many dark satellites do not shine, thus solving the satellite problem, but how and 
wh ich satellites become lumino us is not understood yet. 

Spergel fc Steinhardtl (1200(1 ) suggested another way to solve both problems. They real- 



ized that if dark matter particles would interact with each other with a mean free path of 
the order of the size of galactic cores, the dark matter interactions would efficiently ther- 
malize the system and avoid the formation of both a central cusp and too many satellites. 
The requirement on the mean free path A is roughly A ~ 10 kpc. We can figure out the 
necessary cross section a for dark matter self-interactions by recalling that the mean free 
path is related to the cross section and the number density n, or matter density p = mn, 
through the relationshi p A = l/(na) = ml {pa). T aking a typical p ~ 0.3 GeV/cm'^ gives 



tne reiationsni p A = ilxna] = ml [pa]. l aKmg a typical p ~ u.d (orev/cm" gives 
60 cm^/g. So ISpergel k. Steinhardt ( 2000l ) suggested a new self-interacting dark 



a /m 

matter particle with a /m m. the range 

1 cmVg ^ (j/m ^ 100 cmVg. (32) 

To understand the magnitude of this number, it is useful to compare it with the geometric 
cross section of a proton, which is one of the largest known cross sections for elementary 
particles. We have 

(a/m)proton ~ 1 fmVGeV ~ 0.006 cmVg. (33) 
Thus the desired a/m seems rather big. It is therefore not s urprising that astrophysica l con- 



straints on self-interacting dark matter are rather stringent. iGnedin fc Ostrikeii (|200l[ ) con- 
sidered the evaporation of halos inside cluste rs and set the constraint a/m < 0.3—1 cm^/g. 



Yoshida. Springel. White, fc Tormen ( 2000( ) considered the shape of cluster cores, which is 



rounder for self-interacting dark matter than for standard cold dark matte r, and concluded 
that a /m must be < 10 cm^ /g. This bound was lat er strengthened by Miralda-Escudel 



( 2002h to a/m < 0.02 cm^/g. iMarkevitch et all (j2003[ ) discovered a gas bullet lagging be- 



hind dark matter in the merging galaxy cluster 1E0657-56. They combined Chandra X-ray 
maps of the hot gas in the cluster with weak lensing maps of its mass distribution. From 
estimates of the column mass densities and of the distance between the gas and the dark 
matter, Markevitch et al. were able to set the upper limit a/m < 10 cm^/g with direct 
observations of the dark matter distribution. All these bounds leave little room, if any, to 
self-interacting dark matter a la Spergel & Steinhardt. 
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Figure 8: Relic density of gravitationally-produced WIMPZILLAs as a function of their 
mass Mx- Hi is the Hubble parameter at the end of inflation, Trh is the reheating temper- 
ature, and Mpi ^ 3 x 10^^ GeV is the Planck mass. The dashed and solid lines correspond 
to inflationary models that smoothly end into a radiation or matter dominated epoch, re- 
spectively. The dotted line is a thermal distribution at the temperature indicated. Outside 
the 'thermalization region' W IMPZILLAs cannot reach thermal equilibrium. (Figure from 
Chung. Kolb. fc Riottol . flQQSL ) 



WIMPZILLAs 

Our last example of cold dark matter candidates illustrates a fascinating idea for generating 
matter in the expanding universe: the gravitational creation of matter in an accelerated 
expansion. This mechanism is analogous to the production of Hawking radiation around a 

black hole, and of Unruh radiation in an accelerated reference frame. 

WIMPZILLAs (IChung. Kolb. fc RiottoLll998Lll999l:fKuzmin fc Tkachevl . Eoosh are very 



massive relics from the Big Bang, which can be the dark matter in the universe if their mass 
is ~ 10^'^ GeV. They were produced at the end of inflation through a variety of possible 
mechanisms: gravitationally, during preheating, during reheating, in bubble collisions. It is 
possible that their relic abundance does not depend on their interaction strength but only 
on their mass, giving great freedom in their phenomenology. To be the dark matter today, 
they are assumed to be stable or to have a lifetime of the order of the age of the universe. 
In the latter case, their decay products may give rise to the highest energy cosmic rays, and 
solve the problem of cosmic rays beyond the GZK cutoff. 

Gravitational production of particles is an important phenomenon that is worth describ- 
ing here. Consider a scalar field (particle) X of mass Mx in the expanding universe. Let rj 
be the conformal time and a (77) the time dependence of the expansion scale factor. Assume 
for simplicity that the universe is fiat. The scalar field X can be expanded in spatial Fourier 
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modes as 

(^^^3/2^(^) [^khMe''-' + a\h*Me-^'-'\ . (34) 

Here and a\ are creation and annihilation operators, and hki^r]) are mode functions that 
satisfy (a) the normalization condition h^h'^ — h'^hl = z (a prime indicates a derivative with 
respect to conformal time), and (b) the mode equation 

K{v)+^l{v)hk{v) = 0, (35) 

where 

u;lir^) = e + MW + m-l)-- (36) 



a 



The parameter ,^ is = for a minimally-coupled field and C = | for a conformally-coupled 
field. The mode equation, Eq. ()35j) . is formally the same as the equation of motion of 
a harmonic oscillator with time- varying frequency uJk{f]). For a given positive-frequency 
solution hkirj), the vacuum \0h) of the field X, i.e. the state with no X particles, is defined 
as the state that satisfies a^lOh) = for all k. Since Eq. f^^j) is a second order equation and 
the frequency depends on time, the normalization condition is in general not sufficient to 
specify the positive-frequency modes uniquely, contrary to the case of constant frequency 
coq for which h^{r]) = e~*'^°^/(2ci;o)^''^. Different boundary conditions for the solutions hk{r]) 
define in general different creation and annihilation operators and a\, and thus in general 
different vacua.^ For example, solutions which satisfy the condition of having only positive- 
frequencies in the distant past, 

h{7]) ~ e-^-^fc"^ for r] -oo, (37) 

contain both positive and negative frequencies in the distant future, 

h{r]) ~ ake-"^^"^ + /3fce+^^*^^ for t] +oo. (38) 

Here uj^ = lim^^-too ^^k{v)- a consequence, an initial vacuum state is no longer a vacuum 
state at later times, i.e. particles are created. The number density of particles is given in 
terms of the Bogolubov coefficient in Eq. ()38|) by 



These ideas have been app lied t o gravitational parti c le cre ation at the end of inflation by 
Chung. Kolb. fc Riottol (Il998t ) and lKuzmin fc TkachevI (Il998t ). Particles with masses Mx of 



the order of the Hubble parameter at the end of inflation, Hj 10~^Mpi ^ 10^^ GeV, may 
have been created with a density which today may be comparable to the critical density. 
Figure 8 shows the relic density Qh"^ of these WIMPZILLAs as a function of their mass Mx 
in units of Hj. Curves are shown for inflation models that have a smooth transition to a 
radiation dominated epoch (dashed line) and a matter dominated epoch (solid line). The 

^The precise definiti on of a vacuum in a curved space-time is stil l subjec t to so me ambiguities. We refer 
the interested re ader tolFuUind ((l97ii Il989l) : iBirrell fc Davi^ (|l982l) : Iwaidi (|l994 ) and to the discussion in 



Fchung. Notari. fc Riottol lj2003(l and references therein 
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third curve (dotted line) shows the thermal particle density at temperature T = Hi/2ti. 
Also shown in the figure is the region where WIMPZILLAs are thermal relics. It is clear 
that it is possible for dark matter to be in the form of heavy WIMPZILLAs generated 
gravitationally at the end of infiation. 

3 Neutralino dark matter searches 

The second part of these lectures is an introduction to several methods to detect non- 
baryonic dark matter. We will use the lightest neutralino as our guinea pig, because of 
the variety of techniques that can be employed to detect it, but the discussion is more 
general and can be applied to a generic WIMP. Thus in this second part we assume that 
non-baryonic dark matter is made of WIMPs (in particular, neutralinos) , and we examine 
several observational ways to test our assumption. 

Neutralino dark matter searches are traditionally divided into two main categories: (1) 
direct detection of Galactic dark matter in laboratory experiments, and (2) indirect detection 
of neutralino annihilation products. For the sake of exposition, indirect searches are further 
subdivided into: (2a) searches for high-energy neutrinos from the center of the Sun or of 
the Earth; (2b) searches for anomalous cosmic rays and gamma-rays from galactic halos, 
especially our own; and (2c) searches for neutrinos, gamma-rays, and radio waves from the 
Galactic Center. We now examine each of them in turn. 

3.1 Direct detection 

The idea here is that neutralino dark matter is to be found not only in the halo of our 
galaxy and in our solar system, but also here on Earth and in the room we are in. Thus 
if we could set up a detector that records the passage of dark matter neutralinos, we could 
hope of detecting neutralino dark matter. 

A process that can be used for this purpose is the elastic scattering of neutralinos off 
nuclei. Inelastic scattering could also be used in principle, as could scattering off electrons, 
but the rate of these processes are expected to be (much) smaller. 

Dozens of experiments worldwide, too numerous to be all listed here, are using or plan to 
use elastic scattering to search for neutralino dark matter, or WIMP dark matter in general. 
The small expected detection rate, and the necessity of suppressing any ionizing radiation 
passing through the detector, are reasons to shelter these experiments from cosmic rays, e.g. 
by placing them in mines or underground laboratories. 

Generally, with the notable exception of directional detectors described below, only the 
energy deposited in the detector during the elastic scattering can be measured. This energy 
is of the order of a few keV, for typical neutralino masses and speeds in the galactic halo. 
The kinetic energy of the recoiling nucleus is converted partly into scintillation light or 
ionization energy (giving an electric current) and partly into thermal energy (heating up 
the detector). 

In cryogenic detectors, a simultaneous measurement of both ionization and thermal 
energy allows the discrimination of nuclear recoils from electrons produced in radioactive 
decays or otherwise. This discrimination, however, cannot tell if the nuclear recoil was 
caused by a WIMP or an ambient neutron. The detector, most often a germanium or 
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Figure 9: Annual modulation of the total counting rate (background plus possible dark 
matter signal) in seven years of data with the DAMA-Nal detector. A constant counting 
rate has been subtracted to give the 'residuals.' The significance of the modulation is 6cr 
and its period is 1 year. The inter pretation of the yearly modulation as due to a WIMP 



signal is controversial. (Figure from Bernabei et al. . 20031 ) 



silicon crystal, needs to be cooled at liquid helium temperature so that its low heat capacity 
converts a small deposited energy into a large temperature increase. Only relatively small 
crystals can be currently used in these cryogenic detectors, with relatively low detection 
rates. 

Detection rates can be increased by using bigger detectors operated at room temperature, 
at the expense of giving up a measurement of the thermal energy and loosing discrimination 
power against electrons. The biggest dark matter detector is currently of this type. It is a 
sodium iodide crystal (a scintillator) under the Gran Sasso mountain in Italy, and it belongs 
to the Italian-Chinese collaboration DAMA (short for DArk MAtter). Interestingly, the loss 
of discrimination power and the gain in target mass almost compensate each other, and the 
sensitivity of cryogenic and scintillation detectors is not very different. 



Annual modulation 



Few years ago, the DAMA collaboration reported a possible detection of WI MP dark matter 
(iBelliLllQQTtlBernabei et Z[ll998Lll999[l2000Ll20ol . Their most recent data (|Bernabei et all 
2003| ) span 7 years and show a 6.3a modulation in their total counting rate (signal+background) 
with a period of 1 year and an amplitude of ~ 0.02 events/day per kg of detector and keV 
of visible re coil energy (Figure 9). This kind of y early modulation in a WIMP sig nal was 
predicted by Drukier. Freese. fc Sperge3 ( 19861 ) and lFreese. Frieman. fc Gouldl (|l988l ) on the 
basis that the velocity of the Earth around the Sun adds vectorially to the velocity of the 
Sun in the Galaxy to produce a yearly modulation in the average speed of the WIMPs 
relative to the Earth (the WIMPs are assumed on average at rest in the Galaxy). For an 
observer on the Earth, the WIMP 'wind' arrives at a higher speed when the Earth and the 
Sun move in the same direction and at a lower speed when they move in opposite directions. 
(The two velocities are actually misaligned by ~ 60°; see Figure 10.) The WIMP flux, and 
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the WIMP detection rate, both proportional to the relative speed of the Earth and the 
WIMPs, are similarly modulated. 

While the presence of a yearly modulation in the DAMA data seems now to be estab- 
lished, its interpretation as due to WIMPs is controversial. Firstly, it is not hard to imagine 
that the background itself can undergo seasonal variations with a period of a year. The 
DAMA collaboration has examined many possible sources of background variation, includ- 
ing a yearly modulation of the cosmic ray intensity underground due to winter-summer 
temperature changes in the upper atmosphere. They claim to have found no annual varia- 
tion in the background level that would produce an amplitude as big as the observed one. 
Secondly, the EDELWEISS cryogenic det ector has sensitivity comparable to DAMA's but 
has not recorded any nuclear recoil event ( Benoit et al.1 . 120021 ). And the CDMS-I cryogenic 



detector, also of comparable sensitivity, has detected nuclear recoil events attributed t o 
ambient neutrons in the shallow site where the detector was running ( Akerib et al. . 2003[ ). 



Comparison of these experimental results is however not as straightforward as it may seem, 
because the relationship between detection rates in cryogenic and scintillation detectors de- 
pends, among other things, on the kind of WIMP-nucleus interaction and on details of the 
WIMP velocity distribution in the halo, which are both poorly known. 

This dependence is apparent in the expression for the expected counting rate per recoil 
energy bin and unit detector mass dR/dE. We have 

dR f Nt da , ,q , 

X — X nvf{v,t)d^v, (40) 



dE J Mt dE 



where Nj- and My are the number of target nuclei and the detector mass, respectively, da/dE 
is the WIMP-nucleus differential cross section, and nvf{v,t) is the WIMP flux impinging 
on the detector. Here n denotes the WIMP density, v the WIMP speed, and f{v,t)d^v 
the WIMP velocity distribution. We write Mt/Nj- = M, the nuclear mass, n = p/m, and 
da/dE = (To|-F(g)p/£'max, where ao is the total scattering cross section of a WIMP off a 
fictitious point-like nucleus, |-F(g)p is a nuclear form factor that depends on the momentum 



Figure 10: Sketch illustrating the 
directions of the Sun's and the 
Earth's motions during a year. As 
the Sun moves in the Galaxy (here 
at 232 km/s, 60° out of the plane of 
the Earth's orbit), the Earth moves 
around the Sun (here at 30 km/s). 
The vectorial sum of their velocities 
gives the velocity of the Earth with 
respect to the Galaxy. Assuming 
the WIMPs to be on average at rest 
in the Galaxy, it follows that the av- 
erage speed of the WIMPs relative 
to the Earth is modulated with a 
period of 1 year. 
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transfer q = \/2ME and is normalized as -F(O) = 1, and -Emax = is the maximum 

recoil energy imparted by a WIMP of speed f (/x = mM / {m + M) is the WIMP-nucleus 
reduced mass). Hence 



dR pao|F(g)|2 f f{v,t)^. 



dE 



2m/i Jv>y/ME/2i? 



d'v. (41) 



Notice that one can only measure the product pa^ with this technique. Notice also 
that the event rate at energy E depends on the WIMP velocity distribution at speeds 
V > ^ ME/2p?. This integration limit depends on the nuclear mass, and thus detectors 
with different kinds of nuclei are sensitive to different regions of the WIMP velocity space. 
Moreover, the cross section ctq scales differently for spin-dependent and spin-independent 
WIMP-nucleus interactions. Finally, while there is a consensus on the spin-independent nu- 
clear form factors, spin-dependent form factors are sensi tive to detailed modeling of the pro- 



ton a nd neutron wave functions inside the nucleus (see lJungman. Kamionkowski fc Griest 
19961 and references therein). 



For spin-independent interactions with a nucleus with Z protons and A — Z neutrons, 
one has ^ ^ 

fTo = —\ZGl + {A- Z)^:^ ~ A^^lGPf, (42) 

TT TT 

where and are the scalar four-f ermion couplings of the WIMP with point-like protons 



:ou pnngs c 

and neutrons, respectively (see, e.g, iGondolol Il996t ). The last approximation holds for 



the case ~ G", which is typical of a neutralino. The spin-independent event rate, 
proportional to ao/fi"^, scales with the square of the atomic number A (if we neglect the 
form factor). It is this dependence on A that allows detectors with relatively heavy nuclei 
to reach down to WIMP-proton cross sections typical of weak interactions. 
For spin-dependent interactions, one has instead 

^0 = ——j-{S,)GI + {S,,)G:\ (43) 

where J is the nuclear spin, {Sp) and {Sn) are the expectation values of the spin of the protons 
and neutrons in the nucleus, respectively, and G^ and G" are the axial fqur-fermion coupling s 



ana neutrons m tne nucleus, respectively, and Lr^ and (j„ are tne axial lour-lermion couplmg s 
of the WIMP with point-like protons and neutrons (see Gondolol . 19961: Tovev et al. . 2000[ ). 



There is no increase of the spin-dependent rate with A^, and spin-dependent cross sections 
of the order of weak cross sections are hard to reach with current detector technology. 

Given all these ambiguities in the comparison of cryogenic and scintillator results, it 
is the author's opinion that the important issue if WIMP dark matter has been detected 
is not settled yet. A bigger DAMA detector and an upgraded CDMS detector running in 
the low-background Soudan mine are currently taking data. EDELWEISS is also improving 
their sensitivity, and new experiments, like CRESST-II and ZEPLIN-IV, should start taking 
data shortly. 

Current bounds and future reach 

The sensitivity of some future experiments is shown in Figure 11, togeth er with the current 
best bounds from the cryogenic detectors CDMS-I (|Akerib et al.L and EDELWEISS 
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Figure 11: Current best bounds 
(lines labeled 'CDMS' and 'EDEL- 
WEISS') on the dark matter 
WIMP-proton spin-independent 
cross section clS db function of 
WIMP mass, together with the 
expected reach of a few of the 
upcoming experiments ('CRESST- 
ir, 'CDMS-ir, 'ZEPLIN-IV, 
'XENON'), some theoretical expec- 
tations ('Dirac neutrino', 'Baltz & 
Gondolo', 'CMSSM Ellis et al.'), 
and the first bound on WIMP dark 
matter ('Ahlen et al'). See text for 
details. 



Benoit et al.l. 120021). and the region where DAMA claims evidence for a WIMP signal 
Bernabei et all l2003l ). As it is conventional in comparing results from different experi- 



ments, the figure shows the WIMP-proton spin-independent cross section obtained from 
experimental data using Eq. (j42p under the assumption of a Maxwellian distribution with 
conventional parameters for the WIMP velocity. For an historical perspectiv e, the figure 



also d isplays the first observational bound on WIMP dark matter obtained bv lAhlen et al. 

(mi). 

We also indicate theoretical predictions for a Dirac neutrino with standard model cou- 
plings, a nd for the lightest neu t ralino in two supersymmetric scenarios: minimal supergrav- 
ity as in lEllis. Ferstl. fc Olivel tOO(t (shaded yellow region) and weak-scale MSSM as in 
Baltz fc Gondolol (|200lLl2003l ) (black crosses and magenta squares). Theory models assume 



that the respective dark matter particles fill up the galactic halo. Dirac neutrinos are ex- 
cluded as main constituents of galactic dark matter in the mass range 3 GeV - 3 PeV (the 
bounds continue linearly to the right of the figure). These are the bounds from dark mat- 
ter searches used in Section 2.2 (and Figure 3) to conclude that a yet-undetected particle 
species is needed to provide cold dark matter. For the neutralino, the expected scattering 
cross section varies in a wide range. Even with the most restrictiv e assumptions of the 
Const rained Minimal Supersymmetric Standard Model (CMSSM) in lEllis. Ferstl. &: Olivel 
( ioooh . the cross section at a given neutralino mass can change by an order of magnitude 
when the other supersymmetric parameters are changed (shaded yellow region in Figure 
11). For neutralino masses ~ 400 GeV, the CMSSM model parameters may conspire to 
make the spin- independent cross section arbitrarily small (if such is the case, the total cross 
section will be dominated by spin-dependent terms, which are however much smaller for 
the heavy nuclei in current detectors). The highest CMSSM cross sections are within reach 
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of the upcoming experim ents, although outside the curren t best bounds. Other studies of 
supergravity models (e.g. Feng. Matchev. fc Wilczek . 200Clf ) are less restrictive and result in 
somewhat l arger cross sections. En l arging the parameter space beyond supergravity, as for 
example in Baltz fc Gondolol ( 2001 . 2003[ ). opens up more possibilities for the values of the 
cross section. This happens basically for two main reasons; (i) The supergravity relation 
between the masses of squark and sleptons on one side and Higgs bosons on the other side is 
removed; the Higgs boson can thus be lighter, and the spin-dependent scattering cross sec- 
tion, which varies essentially with the fourth inverse power of the Higgs mass, can be larger; 
(ii) It is no longer required that the electroweak symmetry breaking is achieved radiatively, 
with the consequence that the gaugino and higgsino content of the lightest neutralino can 
be arbitrary; the scattering cross section is then enhanced because mixed neutralinos couple 
to nucleons stronger than p ure gauginos or pure higgsinq s. This explains the larger extent 
of the region covered by the iBaltz fc Gondolol (l200lU2003t ) models in Figure 11. Finally, the 
region in Figure 1 1 mar ked by the magenta squares ind icates the subset of the supersymmet- 
ric models examined bv lBaltz fc Gondolol (I2OOILI2OO3I ) that could be able to quantitatively 
explain the ~ Scr deviation between the measured value and the Standard Model value of 
the magnetic moment of the muon.^ These models have relatively light masses for super- 
symmetric partners, and give neutralino-proton cross sections which are relatively large and 
within the reach of the most ambitious future experiments. 



Directional detection 

We conclude this section by mentioning the very intriguing possibility of WIMP detectors 
that are sensitive to the direction of nuclear recoils (directional detectors). 

The advantages of measuring the recoil direction are multiple: a more powerful back- 
ground discrimination; the detection of the new modulation effects, such as a daily mod- 
ulation in the arrival direction of WIMPs due to the Earth rotation around its axis; and 
the exciting possibility of reconstructing the WIMP velocity distribution in the solar neigh- 
borhood. The latter is possible because of a simple relation between the WIMP velocity 
distributi on f(v,t) and t he nuclear recoil rate differential in both energy E and recoil di- 
rection q ( Gondolol 20021 ) . 



dR nao\F{q)\' 



dEdn 



f{w,q) 



(44) 



where dD. is an infinitesimal solid angle around the direction q, w = ^y ME/2jj?, and 



f{w,q)= / 5{u-q-w)f{v,t)d^v 



(45) 



is the Radon transform of the WIMP velocity distribution. 

A promising development in this direction is the DRIFT detector ( Snowden-Ifft. Martoff. fc Burwell 

2000l ). This detector consists of a negative ion time projection chamber, the gas in the cham- 
ber serving both as WIMP target and as ionization medium for observing the nuclear recoil 

^The size of the deviation has been hard to determine conclusively because of the difficulty of the non- 
perturbative QCD calculations involved, in particular their dependence on the data used as input. There is 
also a story on sign errors in some of the theoretical calculations.... 
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tracks. The direction of the nuclear recoil is obtained from the geometry and timing of the 
image of the recoil track on the chamber end-plates. A 1 m^ prototype has been successfully 
tested, and a 10 m^ detector is under consideration. 

Directional detection is particularly powerful for detecting structure in the dark matter 
velocity space, as discussed in the next section on the Sagittarius stream. 



Sagittarius stream 



Recent observations of the stellar component of the Galactic halo show evidence of a merger 
history that has not yet become well mixed, and corroborate pre vious indications that halos 
form hierarchically. In particular , the Sloan Digital Sky Survev (iNewberg et aL , 20031) and 
the Two Micron All Sk y Survev (IMaiews d. Skrutskie. Weinberg, fc Ostheimed . 2003| ) have 
traced the tidal stream ( Ibata et al.l . 2001 ) of the Sagittarius (Sgr) dwarf spheroidal galaxy. 
The Sagittarius dwarf spheroidal galaxy, of roughly IO^Mq, is a satellite of our own much 
larger Milky Way Galaxy, located insid e the Milky Way, ^ 12 kpc behind the Galactic Center 



and ~12 kpc below the Galactic Plane (jibata et al.Lll997^ . Two streams of matter are being 



tidally pulled away from the main body of the Sgr galaxy and extend outward from it. These 
streams, known as the leading and trailing tidal tails, are made of matter tidally pulled away 
from the Sgr galax y. It appears that the leading tail is showering mat ter down upon the 
solar neighborhood ( Maiewski. Skrutskie. Weinberg, fc Ostheimer . 20031 ). The flow is in the 
general direction orthogonal to the Galactic plane and has a speed of roughly 300 km/s. 
This speed is comparable to that of the relative speed of the Sun and the WIMPs in the 
general dark halo. 

It is natural to expect that dark matter is associated with the detected tidal streams. 
Hence one can hope to detect the stream in direct detection experiments. The detectability 
depends on the density of dark matter in the stream. The mass-to-light ratio M/L in the 
stream is unknown, but is plausibly at least as large as that in the Sgr main body; in fact, 
the M/ L in the stream may be significantly larger because the dark matter on the outskirts 
of the main body would be tidally stripped before the (more centrally located) stars. Various 
determinatio ns of the M/L for the Sgr main body give value s in th e range 25 to 100 (see the 
discussion in Maiewski. Skrutskie. Weinberg. &: Ostheimed. 20031 and refereri c es th erein) . 
Freese. Gondolo. fc Newberd ( 20031 ) and Freese. Gondolo. Newberg. fc Lewid ()2003| ) have 
estimated the density of dark matter in the stream, and find it to be in the range 0.3% 
to 23% of the local (smoothed) dark halo density. This agrees with a previous theoretical 
stu dy on the tidal disruption of sate llite galaxies falling into the halo of our own Milky Way 
bv IStiff. Widrow. fc Frie^^ (I2OOIL These authors found that, with probability of order 
1, the Sun should be situated within a stream of density ~ 4% of the local Galactic halo 
density. 

The additional flux of WIMPs from the stream shows up as a 0.3-23% increase in the 
rate of nuclear recoils at energies below a characteristic energy the highest energy that 
WIMPs in the stream can impart to a target nucleus. Hence, there is a step in the energy 
recoil spectrum; the count rate in the detector is enhanced at low energies, but then returns 
to the normal value (due to Galactic halo WIMPs) at all energies above the critical energy 
Ef.. This feature can be observed as a sharp decrease in the count rate above a characteristic 
energy that depends on the mass of the target nucleus, the mass of the WIMP, and the speed 
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of the stream relative to the detector. Figures 12(a) and (b) show how the recoil spectrum 
dR/dE is modified by the presence of the stream for a sodium iodide detector (like in 
DAMA) and a germanium detector (like for CDMS and EDELWEISS). For the sake of 
illustration, the plots assume a stream density equal to 20% of the local halo density. 

Excitingly, the effect of the stream should be detectable in DAMA, CDMS-II, and other 
upcoming d etectors, and may already be pres ent in the current DAMA data. A detail 
calculation ( Freese. Gondolo. fc Newbere . 20031 ) predicts the presence of stream WIMPs in 



the data with a significance of lOOcr for DAMA and 11 a for CDMS if the stream density is 
20% of the local halo density, and a significance of 24cr for DAMA and Sa for CDMS if the 
stream density is 4% of the local halo density. (These significance figures are however very 
sensit ive to the velocity assumed for the stream, cfr. iFreese. Gondolo. Newberg. fc Lewia . 
200i) 



Directional detection will be a fantastic means of recognizing the presence of a dark 
stream through the Solar system. The recoil distribution due to WIMPs in a stream is 
very much different from the recoil distribution due a Maxwellian velocity distribution. The 
corresponding Radon transforms that appear in Eq. (j44j) are: for a stream of velocity V , 

/(u;,q) ^5(u;-\7-q) (46) 
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Figure 12: Effect of the presence of WIMPs in the Sagittarius leading tidal arm. Count rate 
of 60 GeV WIMPs in (a) an Nal detector such as DAMA and (b) a Ge detector such as 
CDMS and EDELWEISS, as a function of recoil energy. The dotted lines (towards the left) 
indicate the count rate due to Galactic halo WIMPs alone for an isothermal halo. The solid 
and dashed lines indicate the step in the count rate that arises if we include the WIMPs 
in the Sgr stream for Vgtr- = 300 km/s in the direction (/, fe) = (90°, —76°) with a stream 
velocity dispersion of 20 km/sec. The plot assumes that the Sgr stream contributes an 
additional 20% of the local Galactic halo density. The solid and dashed lines are for June 
28 and December 27 respectively, the dates at which the annual modula t ion o f the stream 
is maximized and minimized. (Figure from Freese. Gondolo. fc Newberd . 20031 ) 
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Figure 13: Count rate of 60 GeV WIMPs in a CS2 detector (DRIFT) as a function of recoil 
energy E and direction of the nuclear recoil (qx, Qy, Qz)- Here X points toward the galactic 
center, Y toward the direction of galactic rotation, and Z toward the North Galactic Pole. 
On the left is a density plot of the count rate in the 2-dimensional slice EJqy = —20 keV. 
On the right is the 1- dimensional section through EJqy = —20 keV and qx = 0. The 
horizontal axis represents recoils in the direction of the Galactic center (left) and Galactic 
anticenter (right); the vertical axis represents recoils in the direction of the North Galactic 
Pole (upward) and South Galactic Pole (downward). The gray scale indicates the count 
rate per kilogram of detector per day and per unit cell in the 3-dimensional energy space 
Eq. Lighter regions correspond to higher count rates. The white band in the upper part is 
the location of nuclear recoils due to WIMPs in the Sgr stream. The fuzzy gray cloud at 
the center contains recoils due to WIMPs in the local isothermal Galactic halo. The two 
WIMP popu lations can in principle be eas i ly sep arated, given a sufficiently long exposure. 



epc 

(Figure from lFreese. Gondolo. &: Newberd . l2003l ) 



which is non-zero only on the surface of a sphere in {wq^, wqy, wqz) space; for a Maxwellian 
of bulk velocity V and velocity dispersion a, 



JXw,q) 



V27ra2 
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{w -V ■ qf 
2^2 



(47) 



which is a smooth gaussian distribution. For our Sgr stream, consider a next-generation 
DRIFT detector of 30 m^ (DRIFT-2). The difference in the recoil direction distributions 
of stream and isothermal WIMPs is apparent in Figure 13, where we plot the differential 
detection rate E~'^dR/dEdfl for DRIFT-2 under the assumption of a 20% stream density. 
Aa seen in the figure, a large DRIFT detector will have the capability of clearly identifying 
WIMPs in the Sgr stream. 
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3.2 Indirect detection 

Besides the direct detection of galactic neutralino dark matter in the laboratory, we can 
search for dark matter neutralinos by looking for the products of their annihilation. We 
distinguish three types of searches according to the place where neutralino annihilations 
occur. The first is the case of neutralino annihilation in the Sun or the Earth, which gives 
rise to a signal in high-energy neutrinos; the second is the case of neutralino annihilation 
in the galactic halo, or in the halo of external galaxies, which generates gamma-rays and 
other cosmic rays such as positrons and antiprotons; the third is the case of neutralino 
annihilations around black holes, in particular around the black hole at our Galactic Center. 

All these annihilation signals share the property of being proportional to the square of 
the neutralino density. This follows from the fact that the neutralino is a Majorana fermion, 
i.e. is identical to its antiparticle. Two neutralinos can annihilate to produce standard model 
particles. Simple stoichiometry then tells us that the annihilation rate, being proportional 
to the product of the densities of the initial particles, is proportional to the square of the 
neutralino density. In more detail, we have 



where Fann is the neutralino annihilation rate per unit volume (i.e. the number of neutrali- 
nos that are annihilated per unit volume and unit time), a^nn is the neutralino-neutralino 
annihilation cross section, v is the relative speed of the two annihilating neutralinos, p is the 
neutralino mass density, and m is the neutralino mass. Recall that the annihilation cross 
section (Jann goes as 1/v at small speeds, as required by kinematical arguments, and thus 
the product (Tann^ does not vanish linearly with v (and is not small at the relatively small 
speeds of neutralinos in galactic halos). Notice also that the number of annihilations per 
unit volume and unit time is given by iFann; where the factor of | correctly converts between 
the number of annihilation events and the number of neutralinos that are annihilated (2 per 
annihilation). It is easy to get confused with this factor of |. 

If in the annihilation rate Fann we insert a typical weak interaction cross section and a 
typical value for the average dark matter density in the Universe, the annihilation rate we 
obtain gives undetectably small signals. Indirect detection is possible because dark matter 
is not distributed uniformly in space. Galaxies and clusters of galaxies are overdensities 
in the dark matter field, as is any possible substructure in galactic halos. Furthermore, 
dark matter may be concentrated gravitationally around massive objects, and may even 
get trapped inside planets and stars. The neutralino annihilation rate, proportional to the 
square of the neutralino density, increases substantially in these dark matter concentrations, 
sometimes to the point of giving observable signals. 

3.2.1 High energy neutrinos from the core of the Sun or of the Earth 

Neutralinos floating around the solar system can occasionally collide with nuclei in the 
Sun and in the planets (the Earth, in particular). In these collisions, they may loose 
enough kinetic energy to end up with a speed smaller than the escape speed, thus becoming 
gravitationally trapped. After some time, the trapped neutralinos will sink to the core of 
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Figure 14: Illustration of indirect detection of WIMPs using high-energy neutrinos emitted 
in WIMP annihilations in the core of the Sun or of the Earth. 

the celestial body in which they are captured, and will possibly reach a condition of thermal 
equilibrium (Figure 14). 

Once concentrated in the center, neutralinos annihilate copiously. The annihilation 
rate is maximal when all captured neutralinos annihilate (a condition called equilibrium 
between capture and annihilation). Whether this condition is satisfied depends on the 
relative strength of the annihilation and scatteri ng cross sections, and ultimately on the 
parameters of the particle and halo models. (See I.Tungman. Kamionkowski fc Griest . I199(tI 
and references therein for complete formulas.) 



Figure 15: Principle of operation 
of a water Cherenkov neutrino tele- 
scope. An incoming muon neu- 
trino (here a muon neutrino) is 
converted into a charged lepton (a 
muon) in the material surrounding 
the detector (in the rock at the 
bottom of the sea). The charged 
lepton moves faster than light in 
water and thus Cherenkov light 
(blue cone) is emitted along its 
trajectory. The Cherenkov light 
is collected by the array of pho- 
tomultipliers suspended on strings. 
Cherenkov neutrino telescopes in ice 
work on the same principle. (Back- 
ground figure by Frangois Mon- 
tanet, ANTARES Collaboration; 
tracks and Cherenkov light by the 
present author.) 
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Of the annihilation products produced in the center of the Earth and the Sun only 
the neutrinos make it to the surface; all the other products are absorbed or decay within 
a short distance of production. All three flavors of neutrinos are produced for neutralino 
masses which are currently allowed. Direct production of a neutrino pair is however strongly 
suppressed in neutralino annihilation, due to the Majorana nature of the neutralino. Anni- 
hilation neutrinos are instead produced as secondaries in the decay chains of the primary 
particles produced in the neutralino-neutralino annihilation. As a consequence, the neu- 
trino energy spectrum is a continuum, and the typical energy of neutrinos from neutralino 
annihilations is about a tenth of the neutralino mass. Given the current constraints, this 
means a neutrino energy between few GeVs and few TeVs. 

Neutrinos of this energy can be detected in Cherenkov neutrino telescopes, whose prin- 
ciple of operation is depicted in Figure 15. A charged-current interaction in the material 
surrounding the detector (rock, ice, water) converts the neutrino into its corresponding 
charged lepton, which then radiates Cherenkov light in the detector medium (ice or water). 
Several neutrino telescopes are currently operational (among them the Super-Kamiokande 
detector in Japan and the AMANDA detector at the South Pole), and others are under 
construction or development (IceCube at the South Pole, ANTARES and NESTOR in the 
Mediterranean). Other neutrino telescopes have played a role in dark matter searches in 
the past, such as the 1MB, the Frejus, the MACRO, and the Baksan experiments. 

The current experimental situation for this indirect detection method is summarized in 
Figure 16(a) for neutrinos from WIMPs in the Sun, and Figure 16(b) for neutrinos from 
WIMPs in the Earth. The figures show the current best bounds from the MACRO, Baksan, 
Super-Kamiokande, and AMANDA experiments, as well as the first bound obtained using 
this technique by the 1MB collaboration in 1987. Also shown is the reach of the IceCube 
experiment after an exposure of 10 km^ yr, and the ultimate applicability of this method 
for the Sun, which is set by the emission of high-energy neutrinos in cosmic ray interactions 
with nuclei on the surface of the Sun. 

Expectations from theoretical models in Figure 16 (MSSM with seven weak-scale pa- 
rameters) range from cases which are already excluded by this method to cases which this 
indirect method will be unable to explore. In comparison, direct searches have a different 
coverage of theoretical models. The reach of direct searches is indicated in the Figure 16 
by -|- signs in the left panel and by the shading in the right panel (direct searches exclude 
only some of the theoretical models that are projected onto these regions from the higher- 
dimensional supersymmetric parameter space). There are models that can be explored by 
direct searches and not by indirect searches of high-energy neutrinos from the Sun and 
the Earth. And vice versa, there are theoretical models that cannot be explored by direct 
searches but can be explored by indirect searches of high-energy neutrinos from the Sun or 
the Earth. The latter aspect is illustrated in Figure 17, where models that can be reached 
by indirect searches for neutrinos from neutralinos in the Sun are marked by dots in the 
scattering cross section-mass plane, and compared with the sensitivity of several current 
and future direct search experiments. Several models fall beyond the reach of even the 
most ambitious direct dark matter searches. This shows the (everlasting) complementarity 
between direct and indirect neutralino searches. 
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Figure 16: Indirect searches for neutralino dark matter using high-energy neutrinos from 
(left) the Sun and (right) the Earth. On the vertical axis is the flux of neutrino- 
induced muons that traverse the neutrino telescope, on the horizontal axis is the neutralino 
mass. '1MB87' is the historically first upper limit, 'MACRO,' 'BAKSAN,' 'SuperK,' and 
'AMANDA' are the limits from the corresponding experiments. The regions marked by x 
on the left and by dots on the right are the pred i ctions of supersymmetric models defined 
at the weak scale ( Bergstrom. Edsio. fc Gondolol 19981 : Ahrens et al. . 2002| ). The -|- signs 
on the left and the shaded region on the right indicate the regions where there are models 
that have been excluded by direct dark matter searches. The line labeled 'lOkm^yr' shows 
the maximum reach of such an exposure in IceCube, and the line labeled 'Sun background' 
marks the level of high-energy neutrino emission due to cosmic ray interactions on the sur- 
face of the Sun, whic h is the ultimate applica bility limit of this rn ethod. (Figure on the left 
from lGondoQbnnnI : figure on the right from lAhrens et allbood .) 



3.2.2 Gamma-rays and cosmic rays from neutralino annihilation in galactic 
halos 

We shift now our attention to signals originating in neutralino annihilations which occur in 
the halo of our galaxy or in the halo of external galaxies. 

The annihilation products of importance are those that are either rarely produced in 
astrophysical environments or otherwise have a peculiar characteristic that make them easily 
recognizable. In the first category are rare cosmic rays such as positrons, antiprotons, and 
antideuterons. In the second category are gamma-rays, whose spectrum is expected to 
contain a gamma-ray line at an energy corresponding to the neutralino mass (besides a 
gamma-ray continuum; see Figure 18). The gamma-ray line is produced directly in the 
primary neutralino annihilation into 77 or Z7. Positrons, antiprotons, deuterons, and the 
gamma continuum are generated in the particle cascades that follow the decay of the primary 
annihilation products. Their spectra are therefore broad, with a typical energy which is only 
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a fraction of the neutralino mass, and a shape whose details depend on which annihilation 
channels are dominant. Two neutralinos can in fact annihilate into a variety of primary 
products, depending on their masses and compositions: fermion pairs //, Higgs boson pairs 
HiHj, gauge boson pairs W^W~ , ZZ, etc. 

Detection of the gamma-ray line would be a smoking-gun for neutralino dark matter, 
since no other astrophysical process is known to produce gamma-ray lines in the 10 GeV - 
10 TeV energy range. Good energy resolution is crucial to detect the neutralino gamma-ray 
line. The simulation in Figure 19 shows that the upcoming GLAST detector should have 
an adequate energy resolution. 

For these signals, the dependence of the annihilation rate on the square of the density, 
see Eq. (j48|) . has dramatic consequences. The predicted signals may change by several 
orders of magnitude when the model for the dark matter density is changed, even without 
violating observational bounds on the latter. Truly, these observational limits are not very 
stringent, given the understandable difficulty of measuring the dark matter density. Anyhow 
the problem is currently there, and can be divided roughly into two questions: (1) what is 
the radial dependence of the average dark matter density in a galaxy, especially in our own? 
(2) how much substructure, i.e. clumps and streams, is there in galactic dark halos? 

Dependence on the density profile 

Historically, the density profile of galactic dark halos has been given in terms of empirical 
density profiles whose density is constant in a central region and decreases as at large 
radii. The latter is the main ingredient in obtaining a fiat rotation curve in the outer regions, 
which is a primary evidence for dark matter in galaxies. Central among these functions is 



Figure 17: Complementarity of 
direct and indirect neutralino 
dark matter searches. The 
figure shows that several su- 
persymmetric models that are 
within the (expected) reach of 
a big neutrino telescope of 10 
km^ yr exposure are beyond 
the reach of current and fu- 
ture direct detection experi- 
ments. The vertical axis is 
the product of the neutralino- 
proton spin-independent scat- 
tering cross section dx-p and 
the local neutralino density in 
units of 0.3 GeV/cm^ fccDM- 
The horizontal axis is the neu- 
tralino mas s my. (Figure from 



Duda. Celmini. Condolo. Edajo, Silk - 2002 



Duda et al.l . I2003L ) 
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Figure 18: Two examples of gamma- 
ray spectrum from neutralino anni- 
hilation in the galactic halo. The 
numerical values refer to a specific 
model for the dark halo of our galaxy 
and for observations in the direc- 
tion of the galactic center, but the 
spectral characteristics are general. 
Neutralino annihilations produce a 
gamma-ray line at an energy corre- 
sponding to the mass of the neu- 
tralino, and a gamma-ray continuum 
generated in the particle cascades 
following the primary annihilation. 
The figure illustrates that the shape 
of the continuum spectrum depends 
on the neutralino mass, but notice 
that it also depends on the neu- 
tralino composition. (Figu r e from 
Rergstrom. UHio. Rucklevl. Il 9981.1 



the cored isothermal profile 



Poa" 



J.2 _|_ ^2 



(49) 



where a is called the core radius ( Bahcall fc Soneira . 198Cl( ). This parametrization is so 
simple and so much used that it is sometimes called the 'can onical' densitv profile. Anothe r 
interesting empirical parametrization is the density profile of Persic. Salucci. fc Stell ( 199fi[ ) , 



Pps = Po- 



3(r2 + a2)2 



(50) 
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Figure 19: Simulation of 
a gamma-ray annihilation 
line from the annihilation 
of ~48 GeV neutralinos, su- 
perimposed on a gamma-ray 
background of astrophysical 
origin. The simulation includes 
the finite energy resolution of 
the upcoming GLAST detec- 
tor. (Figure from the GLAST 
Science Brochure.) 
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Figure 20: Dark matter den- 
sity profiles for a galaxy resem- 
bling our own. Models 'BS' 



Bahcall fc Soneiral. 11980 and 'PS 



Persic. Salucci. fc Stel 



199fih 



are 



empirical parametrizations which 
possess a central region with con- 
stant density fcore). Mode ls 'NFW 
( Navarro. Frenk. fc White!. Il99fi ) 



and 'Moore et al' (jMoore et al 
1998[ ) are derived from numerical 
simulations of structure formation 
in the Universe, and in them 
the density in the central region 
increases as a power law of radius 
(cusp). All four models are nor- 
malized to the same total mass and 
virial radius. 
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which provides good fits to rotation curves of hundreds of spiral galaxies (which are not as 
flat as one might think!). 

Numerical simulations of structure formation in the Universe have discovered that pure 
cold dark matter halos do not follow the previous empirical density profiles but instead 
have a universal shape whose parameters depend on the mass (or age) of the system. 
Navarro. Frenk. fc Whit j ( 1996[ ) have found this universal profile to have the form 



^3 

r(r + Ts) 

where the parameter is the radius at which the radial dependence of the density changes 
from to r~^. The empirica l depen dence is then seen as an approximation in the 
transition region. Moore et al. ( 1998| ) suggest instead that the universal profile may be 
steeper at the center, 

PMoore('^) = o/a/ , nq/o • (^2) 

Which of these two profiles better represents the results of numerical simulations is a question 
that must be answered by higher resolution simulations (which seem to be pointing to an 
inner slope 7 that depends on the mass of the system). 

The four profiles mentioned above are plotted in Figure 20 for a galaxy that could be 
our own, with total mass Mvir = 2 x IO^^Mq and virial radius rvir = 428 kpc. 

The essential difference between the empirical and the numerical profiles, for what con- 
cerns neutralino signals, is their behavior at small radii. The empirical profiles have a central 
region with constant density, called a core, while the numerical profiles have a density that 
increases toward the center as a power law, called a cusp. Since the neutralino annihilation 
signals scale as the square of the density, any power law r~'^ with 7 > is bound to give 
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Figure 21: Expected gamma-ray 
flux in the gamma-ray line from 
neutralino annihilation in our 
galactic halo, coming from the 
direction of the Galactic Center 
(without the spike discussed in 
next Section). The photon energy 
is to good approximation equal 
to the neutralino mass. The upper 
set of points is for an NFW m o del 
(iBergstrom. Ullio. fc Buckley . 



1998t ). the lower set for a cored 
isothermal model. The gamma-ray 
fluxes differ by a factor of about 
500. Also shown for approximate 
comparison are the sensitivities of 
current and upcoming gamma-ray 
telescopes to gamma-r ay point 
sources. 



2003) 



(Figure from Gondolol . 



a higher signal from the central region, for a given total mass. Moreover, if 7 > 3/2, the 
annihilation rate in the central region formally diverges, because J p^r'^dr oc J r^'^'^^'^dr. 

It is therefore crucial for our purpose to know which class of proflles, with a core or with a 
cusp, better resembles reality. Unfortunately, constraining a given dark matter proflle using 
the kinematics of the central region is a hard problem, because the dynamics of the central 
regions is usually dominated by the visible matter, and when it is not, like in low surface 
brightness (LSB) spiral galaxies, the central parts are so small that the angular resolution of 
the observations is a major concern. Discrepant data and an apparent lack of universality 
in LSB proflles has generated endless controversies in the astrophysical community. At 
any rate, it must not be forgotten that the proflles mentioned above obtained in numerical 
simulations include cold dark matter only, and astrophysical processes connected for instance 
with the gas and the formation of stars may well modify the density proflle of dark matter 
(theoretical work in this direction is not lacking). 

In the face of this situation, when making predictions for annihilation signals from neu- 
tralino dark matter in the halo, it is prudent for the moment to consider both possibilities, 
core and cusp. 

The considerable differences in indirect neutralino signals between a cored and a cuspy 
density proflle are illustrated in Figure 21 for a gamma-ray signal from neutralino annihi- 
lation from the direction of the Galactic Center (without the spike contribution discussed 
in the next Section). There is a factor of 500 difference between assuming an NFW proflle 
or a cored isothermal proflle in t he theoretical calculation of the ga mma-ray flux in the 
seven-parameter MSSM model in Bergstrom. Ullio. fc Bucklev ( 19981 ). Superposed on the 
plot are the sensitivities to gamma-ray point sources of various gamma-ray telescopes, both 
current and upcoming (as of 2000). The comparison with the theoretical expectations in the 
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figure is not direct, however, because the neutrahno emission may not be point hke in some 
of the telescopes. The sensitivity curves are thus only meant to provide an approximate 
comparison. 

It must be mentioned here that EGRET has detected gamm a-ray emission from the 
direction of the Galacti c Center dMaver-Hasselwander et aL . 1998[ ). However, in a reanalysis 



of the EGRET data bv lHooper fc Dinguj (1200^ the EGRET signal seems to orig i nate f rom 



a source which is displaced with respect to the Galactic Center. Hooper fc DingusI ( 20021 ) use 



their reanalysis to place constraints on the gamma-ray flux from neutralino annihilations, 
coming to an upper limit on the gamma-ray flux above 1 GeV which is about a factor of 
2 higher than the theoretical predict i ons fo r neutralino masses ~ 50 GeV, assuming an 
NEW profile. See iHooper fc DingusI (|2002h for details. 



Effect of halo substructure 

The same numerical simulations that predict a cuspy dark halo profile also predict the 
existence of numerous dark clumps in galactic halos. These clumps are a natural outcome 
of the hierarchal formation of structure in the Universe. Small structures form first, and 
larger structures, like galaxies, grow by attracting and swallowing smaller structures. This 
process continues to the present day. Clumps that fall into a galaxy are pulled apart by 
tidal interactions: the material pulled out forms tidal streams that crisscross the galactic 
halo. The central parts of some of the clumps may survive for a long time, and become a 
galactic satellite. 

The overall picture of hierarchical structure formation has found confirmation in a va- 
riety of context, from observations of galaxy clusters and of merging galaxies to the halo 
substructure detected in our own galaxy (see the discussion on the Sagittarius stream in 
Section 3.1 above, for example). A numerical discrepancy should however be mentioned. A 
counting of visible satellites of our own galaxy gives a number of luminous satellites that is 
much smaller than the expected number of dark satellites. A resolution to this discrepancy 
may be that only a small fraction of dark satellites becomes luminous. It is not clear why 
this should happen, but on the other hand we do not fully understand how galaxies become 
luminous in the first place. Thus for what concerns signals from neutralino annihilation, 
it makes sense to examine the effect of adding substructure, i.e. clumps and streams, to 
galactic dark halos. 

Substructure in the halo tends to increase the annihilation signals because of the de- 
pendence of the annihilation rate on the square of the dark matter density (see Eq. (|48|) ) . 
The enhancement factor is linearly proportional to the density enhancement 6 = p' / p. To 
understand why the dependence is linear instead of quadratic in 6, consider a box of volume 
V containing a total mass M. The density in the box is p = M/V and the annihilation rate 
integrated over the whole box is 

-Kami = Tannl/ = 5 7^- (53) 

V 

Now let all the mass be concentrated equally into N small boxes, each of volume V, so 
that each box contains a mass M' = M/N and thus has a density p' = M' /V. The density 
enhancement is then 6 = V/NV. The annihilation rate from the whole box is now given 
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Figure 22: Expected isotropic 
gamma-ray background from neu- 
tralino annihilations in the early 
Universe. Dotted lines are the 
signal from neutralino annihila- 
tions; solid lines are the sum of the 
neutralino signal and a gamma-ray 
background of astrophysical origin. 
Two neutralino models are shown, 
together with the current EGRET 
measurements of the isotropic 
extragalactic background. The as- 
sumptions used for the dark matter 
profile are indicated. (Figure fro m 



Ullio. Bergstrom. Edsio. fc Lacev 



by a sum over the TV small boxes as 

Knn = NT^^^V = N —-ryr = —-^^RTJ = ^^^'^^^ (^4) 

m"^ V m'^ V JSV' 

One power of the density p in Fann is compensated by a decrease in the volume where 
annihilations occur. Hence the signal enhancement is linear in the density increase. 

An interesting consequenc e of th e an nihilation signal enhancement has been explored 
by Bergstrom. Edsjo. fc UlliJ ( 2001 ) and Ullio. Bergstrom. Edsjo. fc Lacev ( 2002[ ). These 



authors have found that the density enhancements produced during the formation of the 
large scale structure in the Universe may lead to a substantial increase in the isotropic 
gamma-ray background from neutralino annihilations in the early Universe. Moreover, they 
found that this increase in the gamma-ray signal is not very sensitive to details of the galactic 
density profile. Expected gamma-ray spectra may in some models be close to the measured 
gamma-ray background, as illustrated in Figure 22. The gamma-ray spectra include both a 
continuum part and gamma-ray lines (two for each neutralino case: one for XX ~^ 77? the 
other for xX ^l)- The gamma-ray lines are asymmetrically broadened because photons 
emitted at earlier times have a larger redshift. The gamma- ray background due to neutralino 
annihilations should be searched for at high galactic latitudes, where the galactic emission 
is expected to be minimal. Detection of the line features depicted in Figure 22 would not 
require an energy resolution much better than the present one. 

Another exciting application of a signal enhancement due to clumps in the galactic halo 
is the boost of the positron signal from neutralino annihilation up to the level of the excess 
of cosmic ray positrons observed by the HEAT collaboration. The HEAT collaboration 
flew two different detectors on balloons, and claims to have detected a ratio of positron to 
electron fluxes above ~7 GeV that is higher than the flux expected in current state-of-the- 
art models of cosmic ray production and propagation in the galaxy. These models aim at 
explaining all correlated signals in gamma-rays, radio waves, protons, electrons, positrons. 
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heavy nuclei, etc. that are produced by cosmic rays in our galaxy. 

One possible explanation is that the positron excess is due to the positrons generated 
in neutralino annihilation in the galactic halo. If the neutralinos are produced thermally 
in the early universe, which is the most common assumption, the annihilation cross sec- 
tion (Tann IS forced to bc small by the requirement that the neutralino relic density is large 
enough for neutralinos to be the dark matter. Using the average value of the local dark 
matter density, of the order of 0.3 GeV/cm^, then leads to a positron signal which is more 
than an order of magnitude smaller than the excess measured by HEAT. Increasing the 
annihilation cross section (Xann does not make the signal higher, because the density p de- 
creases inv ersely with (Tann, and hence t he annihilation signal, being proportional to CannP^, 
decreases. Kane. Wang, fc Wellsl » suggested that ti>e ueutralinos may not have been 
produced thermally in the early Universe, and were thus able to decouple the annihilation 
rate in t he halo from the constraint cora i iig ou t of the relic density requirement. Alter- 
natively, iBaltz. Edsio. Freese. fc Gondolol tm± have suggested that substructure in the 
galactic halo may provide the necessary boost factor to the positron signal. 

Enhancing the positron signal throu gh clumps also enhances other armihilat ion signals, 
such as antiprotons and gamma-rays. iBaltz. Edsio. Freese. fc Gondolol tm± have per- 
formed a detailed analysis of these enhancements, and have concluded that it is possible to 
explain the HEAT positron excess with boost factors as small as 30, but typically higher. 
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Figure 23: Two examples of neutralino models that provide a good fit to the excess of cosmic 
ray positrons observed by the HEAT collaboration. The two sets of data points (open and 
filled squares) are derived from two different instruments flown in 1994-95 and 2000. The 
lines represent: (i) the best expectation we have from models of cosmic ray propagation in 
the galaxy ('bkg. only fit'), which underestimate the data points above ~7 GeV; (ii) the effect 
of adding positrons from neutralino annihilations (lines 'SUSY component', 'SUSY-|-bkg. 
fit', and 'bkg. component', the latter being the resulting background component when the 
data are fitted to the sum of backgro und and neutralino contributions). (Figures from 
Baltz. Edsi5. Freese. fc Gondolol 12002 ) . 
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Figure 24: Sensitivity of upcoming gamma-ray telescopes to neutralino models that can 
explain the HEAT positron excess with neutralino clumps in the galactic halo. Model points 
are indicated by crosses; circles denote those models that in addition can also account for 
the measured deviation in the muon magnetic moment. The upper set of sensitivity curves 
corresponds to the high latitude gamma-ray line flux (scale on the left); the lower set of 
curves to the direction toward the galactic cente r (scale on the right; no steep spike arou nd 
the central black hole is assumed). (Figure from Baltz. Edsio. Freese. fc Gondolol 2002i ) 

without obtaining too many antiprotons or gamma-rays. Figure 23 shows two examples of 
neutralino models that provide a good fit to the positron excess. On the left, the neutralino 
has mass = 340 GeV and is an almost pure gaugino (gaugino fraction Zg = 0.98); the 
boost factor is 95 and other parameters are listed in the figure. On the right, the neutralino 
is mixed (gaugino fraction Zg = 0.70) with a mass of 238 GeV; the boost factor is 116.7. 
The degree of freedom is quite good for both fits, 1.34 and 1.38 respectively. 

The ultimate test of the explanation of the positron excess by means of neutralino 
clumps will be the detection of a signal in gamma-rays. Gamma-ray production would 
in fact be enhanced by the same me chanism that would enhance positron production. 
Baltz. Edsio. Freese. &: Gondolol ( 20021 ) have found that almost all neutralino models that 
can explain the positron excess are within the sensitivity reach of upcoming gamma-ray tele- 
scopes (see Figure 24). The realistic possibility of confirming (or disproving) the neutralino 
origin of the positron excess is fascinating. 



3.2.3 Signals from neutralino annihilation at the Galactic Center 

The last indirect neutralino signals we consider are neutrinos, gamma-rays and radio waves 
from a possible dark matter concentration around the black hole at the Galactic Center. 

Evidence for the presence of a black hole at the center of our galaxy comes from studies 
of the motion of stars in orbit around the center. The speeds of these stars decrease from 
the center as the inverse square root of the radius, which is the primary indication for 
the existence of a point mass at the center. The mass of the central object is measured 
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Figure 25: Density profiles of spikes that form adiabatically around the black hole at the 
center of our Galaxy. The position of the Sun i s indicated by a cross. Four m odels for the 
halo profile are shown: t wo with cores ('PS' bv lPersic. Sa ucci. fc Stell (Il996h and 'can' by 
Bahcall fc Soneiral (Il98nl)) and tw o with cusps f'NFW bv lNavarro. Frenk. fc Whitel (|l996l ) 
and 'Mfc' bv iMoore et ah ( 1998[ )). The spikes form within the radius of influence of the 
black hole, rinfl ~ 1 pc. In the 'annihilation plateau' neutra lino annihilations hav e been so 



rapid as to deplete the number of neutralinos. (Figure from Buckley et al. . 200 ll ) 



to be ^ 4 X 10^ solar masses, which are contained within a sphere of less than ^ 0.05 
DC flEckart fc Genzel \l99^ iGhez. Klein. Morris, fc Beckhnl . Il998l: iGhez et al.L liool . No 

stellar or gas system can be so dense, indicating that the central object is most probably a 
black hole. The position of the black hole happens to coincide with the position of a strong 
radio source called Sagittarius A*, which is thus identified with the central black hole. 

The radio emission from Sgr A* is easily explained by thermal emission of hot matter 
falling into the black hole. However, contrary to many of the similar black holes observed at 
the center of external galaxies, our galactic black hole does not emit intensely in the X-ray 
band, and it is controversial if it emits gamma-rays. Models for such 'quiet' black holes do 
exist, however, such as those involving advection-dominated accretion flows (ADAFs). 

Further evidence for a black hole at the center of the Milky Way comes from th e 2001 
observation of a X-ray and infrared flares from the Galactic Center ( Baganoff et al.l . |2001 



Porauet et all booi iG^zel et all l2003bf) . The flare time scale and intensity can nicely be 



expla ined if the flare is produced near a black hole (see, e.g., Aschenbach. Grosso. Porquet. fc Predehll . 



200J. 



Dark matter may be driven near the black h ole gravitationally, and may form a dense 



concentration around it ( Gondolo fc Silk . 1999f ). We will call this concentration a spike, 
so as to distinguish it from the dark matter cusps of Section 3.2.2. The formation of a 
spike is gravitational phenomenon similar to but less efficient than accretion of matter, in 
that the latter involves dissipation of energy and angular momentum and can thus produce 
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concentrations which are smaller and denser. 

How strong is the dark matter concentration around Sgr A*, or around a generic black 
hole? This is still a matter of investigation . The simplest case is that of adiabatic compres- 
sion, and was analyzed by Gondolo fc Silk ( 1999| ). 

Adiabatic compression of an initial dark matter distribution produces two kinds of spikes. 
If the initial distribution before the black hole forms is thermal, the spike is shallow, with 
density profile p oc r~^/^. If the initial phase-space distribution is a power-law in energy, the 
spike is steep, i.e. p oc with 7 > 2. The physical reason for the higher concentration in 
the second case is the presence of many dark matter orbits with low speed, which are more 
easily driven into bound orbits when the black hole forms. 

These two kinds of spike are illustrated in Figure 25 f or the same four halo rn odels 
discussed in Section 3. 2 .2. M odels with a core, like those by Bahcall fc Soneira ( IQSOl ) and 
Persic. Salucci. fc Stell ( 1996[ ). give rise to a shallow spike around the c entra l black hole, 
while models with a cusp, like those of Navarro. Frenk. fc Whit j f 1996t ) and Moore et al. 
( 1998[ ). produce steep spikes. In the very inner regions, the density may become so high 



that neutralino-neutralino annihilations may have had the time to deplete the number of 
neutralinos and an 'annihilation plateau' is formed. The typical radius of a spike around a 
black hole is determined by the radius of influence of the black hole, rinfi ~ GM/cr^, which 
is the radius at which the gravitational potential energy becomes comparable to the typical 
kinetic energy in the dark matter gas (M is the black hole mass and (T„ is the gas velocity 
dispersion). For the black hole at the Galactic Center, rjnfi ~ 1 pc. 

Neutralino annihilation is enhanced in the spike, because of the dependence of the an- 
nihilation rate on square of the density. The Galactic Center then becomes a source of 
neutrinos, gamma-rays, radio waves, etc. from neutralino annihilation (Figure 26). The 
intensity of the emission depends on the steepness of the spike. If the spike is shallow, neu- 
tralino annihilation is generally undetectable. On the contrary, a steep spike at the Galactic 
Center produces interesting signals. 

For example, if, disregarding Hooper & Dingus's reanalysis mentioned above, we at- 
tribute the EGRET gamma-ray emission from the Galactic Center to neutralino annihila- 
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Figure 26: Artistic conception of 
emission from a dark matter spike 
around a black hole at the galactic 
center. Neutralino annihilation in 
the spike produces intense fluxes of 
neutrinos, gamma-rays, radio waves, 
etc. Some of these signals may be 
detectable. 
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tion in a spike born out of an NFW profile, we obtain a high-energy neutrino flux that is 
either excluded or mostly detectable in a km^ neutrino telescope (Figure 27). The flux of 
neutrino-induced muons above 25 GeV would be detectable over the atmospheric neutrino 
background for neutralino masses between ~100 GeV and ~2 TeV, while heavier neutrali- 
nos would alrea dy be excluded from the current limit on the neutrino emission from the 



Galactic Center (jHabig et all l200lh 



As a second and more dramatic example ( Gondolol 2000bl ). electrons and positrons from 
neutralino annihilation would emit synchrotron radiation as they spiral in the magnetic 
field that plausibly exists around the central black hole. While this synchrotron radiation 
is innocuous for a shallow adiabatic spike, it may exceed the observed radio emission by 
several orders of magnitude if the spike is steep. The radio synchrotron emission at 408 
MHz is shown in Figure 28 for an adiabatic spike born out of an NFW profile, under two 
assumptions for the radial dependence of the magnetic field (a constant field of 1 niG and 
a field in equipartition with the infalling gas). In both cases, all dark matter neutralinos in 
the seven-parameter MSSM models considered are strongly excluded. 

The examples above have assumed that the spike formed adiabatically and maintained 
its shape till the present time. This may not be the case. 

For example, if the central black hole formed through the merging of two black holes 
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Figure 27: Predicted neutrino signal from neutralino dark matter annihilation in a steep 
adiabatic spike at the Galactic Center expressed as number of neutrino-induced muons per 
km^-yr in a neutrino telescope. The spike corresponds to an NFW profile. Each dot in the 
figure corresponds to a point in a seven-parameter weak-scale MSSM, and is normalized 
so that the gamma-ray flux from the spike coincides with the gamma-ray signal from the 
Galactic Center observed by EGRET. Models with heavy neutralinos are excluded by the 
current limit on the neutrino emission from the Galactic Center; models with neutralinos as 
light as ^100 GeV co uld be detected above the atmospheric neutrino background. (Figure 
from iGondoloL l2nn2bL ) 
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Figure 28: Synchrotron emission 
at 408 MHz expected from neu- 
tralino dark matter annihilation in a 
steep adiabatic spike at the Galac- 
tic Center. The spike corresponds 
to an NFW profile, and the syn- 
chrotron radiation is emitted by 
electrons and positrons produced 
in neutralino annihilation. Up- 
per panel: constant magnetic field; 
lower panel: equipartition magnetic 
field. All dark matter neutralinos 
in the seven-parameter weak-scale 
MSSM considered are excluded by 
several o rders of rnagnitu de. (Fig- 
ure from iGondoloL bnOObl ) 
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of similar mass, Merritt. Milosavlievic. Verde, fc Jimenez! ( 20021 ) have shown that the spike 
would become shallow at the end of the merging, because dark matter particles would be 
kicked out of the black hole region via a gravitational sling-shot effect. The final shallow 
spike would not give dramatic signals from neutralino annihilation at the Galactic Center. 
In a realistic scenario of this type, the two merging black holes of similar mass would be 
accompanied by their host galaxies whose mass would also be similar. The merging would 
then constitute what is known as a major merging. A major merging is capable of destroying 
the Galactic disk, and so must not have occurred after the disk formed about 10^° years 
ago. Thus, for this scenario to work, black holes of millions of solar masses should already 
have been in place at very early times. Is this possible? While hard to explain theoretically, 
supermassive black holes have indeed been observed in very distant quasars, at redshift 
2; > 6, so the scenario may be plausible. 

In these considerations, it must be kept in mind that there is a stellar spike around the 
black hole at the Galactic Center. T he steepness o f this stellar spike is however not very 
well know. With large uncertainties, Genzel et al. ( 2003| ) estimate the slope of the stellar 
spike to be 7stars ~ 1.3-1.4. This means that the current stellar spike is probably shallow. 
We may think that the stellar spike is our best proxy for the dark matter spike. If so, also 
the dark matter spike would also be shallow, and thus inconsequential for neutralino signals. 
However, the dark matter and stellar spikes follow very different evolution histories, because 
contrary to the dark matter, binary collisions of stars and coalescence of two stars into one 
at collisions effectively relax the stellar system to a shallower spike. 

The final word has not yet been said regarding the distribution of dark matter at the 
Galactic Center, or around black holes and other compact objects in general. This is one 
of the exciting points of contact between the study of dark matter and the study of the 
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formation and evolution of galactic nuclei. 



4 Conclusions 

Current cosmological data imply the existence of non-baryonic dark matter. We have dis- 
cussed some of the most popular candidates and shown that none of the candidates known 
to exist, i.e. the active neutrinos, can be non-baryonic cold dark matter. Hence to explain 
the nature of cold dark matter we need to invoke hypothetical particles that have not been 
detected yet. Some of these hypothetical particles have been suggested for reasons different 
than the dark matter problem (such as sterile neutrinos, neutralinos, and axions), some 
others have been proposed mainly as a solution to the cold dark matter problem (e.g., 
self-interacting dark matter, WIMPZILLAs, etc.). Although most studies focus on the first 
category of candidates, especially neutralinos and axions, we should keep an open mind. 

To illustrate how we can find out if dark matter is made of elementary particles, we 
have used neutralino dark matter as our guinea pig to survey several methods to search 
for non-baryonic dark matter. These methods range from a direct detection of dark matter 
particles in the laboratory to indirect observation of their annihilation products produced 
in the core of the Sun or of the Earth and in galactic halos, including our own. Direct 
searches may have found a signal from WIMPs (the annual modulation), but this claim 
is highly controversial at the moment. Future direct searches have great promise, and 
might even be able to explore the local velocity distribution of WIMPs. These searches 
are complemented by indirect searches for high-energy neutrinos from the core of the Sun 
or of the Earth. Indirect searches using gamma-rays and cosmic rays from annihilations 
in the galactic halo are subject to uncertainties related to the detailed structure of the 
dark matter halo. Even more so are predictions for dark matter signals from the Galactic 
Center. Despite this, some anomalies in cosmic ray fluxes, namely a positron excess, may 
be explained by neutralino annihilation, and future gamma-ray observations may discover 
a gamma-ray line from neutralino annihilation in our galactic halo. 

All of the examples we have presented are without doubt simple, elegant, and compelling 
explanations for the nature of non-baryonic dark matter. As we ponder on which one of 
them is realized in Nature, we must remember the words of astrophysicist Thomas Gold (as 
quoted by Rocky Kolb): "For every complex natural phenomenon there is a simple, elegant, 
compelling, wrong explanation." 

References 

Abazajian, K., Fuller, G. M., & Patel, M. 2001. "Sterile neutrino hot, warm, and cold dark 
matter," Phys. Rev. D64, 023501 

Abrams, D., et al. [CDMS Collaboration] 2002. "Exclusion hmits on the WIMP nucleon 
cross-section from the cryogenic dark matter search," Phys. Rev. D66, 122003 

Ahlen, S. P., Avignone, F. T., Brodzinski, R. L., Drukier, A. K., Gelmini, G. & Spergel, 
D. N. 1987. "Limits On Cold Dark Matter Candidates From An Ultralow Background 
Germanium Spectrometer," Phys. Lett. B195, 603 



44 



Ahmad, Q. R., et al. [SNO Collaboration] 2002. "Direct evidence for neutrino flavor trans- 
formation from neutral-current interactions in the Sudbury Neutrino Observatory," Phys. 
Rev. Lett. 89, 011301 

Ahrens, J., et al. [AMANDA Collaboration] 2002. "Limits to the muon flux from WIMP 
annihilation in the center of the earth with the AMANDA detector," Phys. Rev. D66, 
032006 

Akerib, D. S., et al. [CDMS Collaboration] 2003. "New results from the Cryogenic Dark 
Matter Search experiment," Phys. Rev. D68, 082002 

Allanach, B. C, Kraml, S. & Porod, W. 2003. "Theoretical uncertainties in sparticle mass 
predictions from computational tools," JHEP 0303, 016 

Aschenbach, B., Grosso, N., Porquet, D., & Predehl, P. 2004. "X-ray flares reveal mass and 
angular momentum of the Galactic Center black hole." Preprint astro-ph/0401589 

Asztalos, S. J., et al. 2004. "An improved RF cavity search for halo axions," Phys. Rev. 
D69, 011101 

Baganoff, F. K. et al. 2001. "Rapid X-ray flaring from the direction of the supermassive 
black hole at the Galactic Centre," Nature 413, 45 

Bahcall, J. N. & Soneira, R. M. 1980. "The universe at faint magnitudes. I - Models for the 
galaxy and the predicted star counts," Astrophys. J. Suppl. 44, 73 

Baltz, E. A., Edsjo, J., Freese, K., & Gondolo, P. 2002. "Cosmic ray positron excess and 
neutralino dark matter," Phys. Rev. D65, 063511 

Baltz, E. A. & Gondolo, P. 2001. "Implications of Muon Anomalous Magnetic Moment for 
Supersymmetric Dark Matter," Phys. Rev. Lett. 86, 5004 

Baltz, E. A. & Gondolo, P. 2003. "Improved constraints on supersymmetric dark matter 
from muon g-2," Phys. Rev. D67, 063503 

Belli, P. 1997. Talk at TAUP 97, LNGS, Italy [published in R. Bernabei et al., Nucl. Phys. 
B (Proc. Suppl.) 70 (1999) 79]. 

Benoit, A., et al. [EDELWEISS Collaboration] 2002. "Improved exclusion limits from the 
EDELWEISS WIMP search," Phys. Lett. B545, 43 

Bergstrom, L. 2000. "Non-baryonic dark matter: observational evidence and detection meth- 
ods," Rept. Prog. Phys. 63, 793 

Bergstrom, L., Edsjo, J., & Gondolo, P. 1998. "Indirect detection of dark matter in km-size 
neutrino telescopes," Phys. Rev. D58, 103519 

Bergstrom, L., Edsjo, J., & Ullio, P. 2001. "Spectral gamma-ray signatures of cosmological 
dark matter annihilations," Phys. Rev. Lett. 87, 251301 



45 



Bergstrom, L. & Gondolo, P. 1996. "Limits on direct detection of neutralino dark matter 
from b — > S7 decays," Astropart. Phys. 5, 263 

Bergstrom, L., UUio, P., & Buckley, J. H. 1998. "Observability of gamma rays from dark 
matter neutralino annihilations in the Milky Way halo," Astropart. Phys. 9, 137 

Bernabei, R., et al. 1998. "Searching For Wimps By The Annual Modulation Signature," 
Phys. Lett. B424, 195 

Bernabei, R., et al. 1999. "On a further search for a yearly modulation of the rate in particle 
dark matter direct search," Phys. Lett. B450, 448 

Bernabei, R., et al. 2000. "Search for WIMP annual modulation signature: Results from 
DAMA/NaI-3 and DAMA/NaI-4 and the global combined analysis," Phys. Lett. B480, 
23 

Bernabei, R., et al. 2003. "Dark matter search," Riv. Nuovo Cim. 26, 1 

Binetruy, P., Girardi, G., & Salati, P. 1984. "Constraints On A System Of Two Neutral 
Fermions From Cosmology," Nucl. Phys. B237, 285 

Birrell, N. D. & Davies, P. C. W., 1982. "Quantum Fields in Curved Space" (Cambridge: 
Cambridge University Press) 

Buckley, J., et al. 2001. "Gamma-ray summary report," in Proc. of the APS/DPF/DPB 
Summer Study on the Future of Particle Physics (Snowmass 2001) , ed. N. Graf, eConf 
C010630, P407 astro-ph/0201160j 

Cembranos, J. A. R., Dobado, A., & Maroto, A. L. 2003. "Brane-world dark matter," Phys. 
Rev. Lett. 90, 241301 

Cheng, H. C, Feng, J. L., & Matchev, K. T. 2002. "Kaluza-Klein dark matter," Phys. Rev. 
Lett. 89, 211301 

Chung, D. J. H., Kolb, E. W.., & Riotto, A. 1988. "Nonthermal supermassive dark matter," 
Phys. Rev. Lett. 81, 4048 

Chung, D. J. H., Kolb, E. W.., & Riotto, A. 1999. "Superheavy dark matter," Phys. Rev. 
D59, 023501 

Chung, D. J. H., Notari, A., & Riotto, A. 2003. "Minimal theoretical uncertainties in infla- 
tionary predictions," JCAP 10, 012 

Dine, D., Fischler, W., & Srednicki, M. 1981. "A Simple Solution To The Strong CP Problem 
With A Harmless Axion," Phys. Lett. B104, 199 

Dolgov, A. D. 2002. "Neutrinos in Cosmology," Phys. Rep. 370, 333 

Drukier, A. K., Freese, K., & Spergel, D. N. 1986. "Detecting Cold Dark Matter Candi- 
dates," Phys. Rev. D33, 3495 



46 



Duda, G., Gelmini, G., Gondolo, P., Edsjo, J., & Silk, J. 2003. "Indirect detection of a 
subdominant density component of cold dark matter," Phys. Rev. D67, 023505 

Edsjo, J. & Gondolo, P. 1997. "Neutralino relic density including coannihilations," Phys. 
Rev. D56, 1879 

Edsjo, J., Schelke, M., Ullio, P., & Gondolo, P. 2003. "Accurate relic densities with neu- 
tralino, chargino and sfermion coannihilations in mSUGRA," JCAP 0304, 001 

Eckart, A., & Genzel, R. 1997. "Stellar proper motions in the central 0.1 pc of the Galaxy," 
MNRAS 284, 576 

Ellis, J. R., Ferstl, A., & Olive, K. A. 2000. "Re-evaluation of the elastic scattering of 
supersymmetric dark matter," Phys. Lett. B481, 304 

Ellis, J. R., Hagelin, J. S., Nanopoulos, D. V., Olive, K. A., & Srednicki, M. 1984. "Super- 
symmetric Relics From The Big Bang," Nucl. Phys. B238, 453 

Feng, J. L., Matchev, K. T., & Wilczek, F. 2000. "Neutralino dark matter in focus point 
super symmetry," Phys. Lett. B482, 388 

Freese, K., Frieman, J. A., & Gould, A. 1988. "Signal Modulation In Cold Dark Matter 
Detection," Phys. Rev. D37, 3388 

Freese, K., Gondolo, P., & Newberg, H. 2003. "Detectability of weakly interacting massive 
particles in the Sagittarius dwarf tidal stream." Preprint astro-ph/0309279 

Freese, K., Gondolo, P., Newberg, H., & Lewis, M. 2003. "The Effects of the Sagittarius 
Dwarf Tidal Stream on Dark Matter Detectors," Phys. Rev. Lett., to appear [Preprint 
astro-ph/ 0310334| 

Fukuda, Y., et al. [Super-Kamiokande Collaboration] 1998. "Evidence for oscillation of 
atmospheric neutrinos," Phys. Rev. Lett. 81, 1562 

Fukugita, M., Hogan, C. J., & Peebles, P. J. E. 1998. "The Cosmic Baryon Budget," Astro- 
phys. J. 503, 518 

Fulling, S. A. 1979. "Remarks On Positive Frequency And Hamiltonians In Expanding 
Universes," Gen. Rel. Grav. 10, 807 

Fulling, S. A., 1989. "Aspects of Quantum Field Theory in Curved Spacetime" (Cambridge: 
Cambridge University Press) 

Genzel, R. et al. 2003. "The Stellar Cusp Around the Supermassive Black Hole in the 
Galactic Center," Astrophys. J. 594, 812 

Genzel, R., Schodel, R., Ott, T., Eckart, A., Alexander, T., Lacombe, F., Rouan, D., & 
Aschenbach, B. 2003b. "Near-infrared flares from accreting gas around the supermassive 
black hole at the Galactic Centre," Nature 425, 934 



47 



Ghez, A. M., Klein, B. C, Morris, M., & Becklin, E. E. 1998. "High Proper-Motion Stars 
in the Vicinity of Sagittarius A*: Evidence for a Supermassive Black Hole at the Center 
of Our Galaxy," Astrophys. J. 509, 678 

Ghez, A. M., Salim, S., Hornstein, S. D., Tanner, A., Morris, M., Becklin, E. E., & 
Duchene, G. 2003. "Stellar Orbits Around the Galactic Center Black Hole." Preprint 
astro-ph/0306130 

Gnedin, O. Y. & Ostriker, J. P. 2001. "Limits on CoUisional Dark Matter from Elliptical 
Galaxies in Clusters," Astrophys. J. 561, 61 

Goldberg, H. 1983. "Constraint On The Photino Mass From Cosmology," Phys. Rev. Lett. 
50, 1419 

Gondolo, P. 1996. "Phenomenological introduction to direct dark matter detection," in 
XXXI Rencontres de Moriond: Dark Matter in Cosmology, Quantum Measurements, Ex- 
perimental Gravitation, Les Arcs, France hep-ph/9605290| . 

Gondolo, P. 2000. "Indirect detection of dark matter." Plenary talk at Neutrino 2000, 
Sudbury, Canada. 

Gondolo, P. 2000b. "Either neutralino dark matter or cuspy dark halos," Phys. Lett. B494, 
181 

Gondolo, P. 2002. "Recoil momentum spectrum in directional dark matter detectors," Phys. 
Rev. D66, 103513 

Gondolo, P. 2002b. "Neutrinos from dark matter." Talk at the NSF Meeting on Neutrinos 
and Subterranean Science, Washington, D.C. 

Gondolo, P., Edsjo, J., Bergstr5m, L., UUio, P., & Baltz, E. A. 2000. "DarkSUSY: A 
numerical package for dark matter calculations in the MSSM." Preprint astro-ph/0012234 

Gondolo, P., Edsjo, J., Ullio, P., Bergstrdm, L., Schelke, M., & Baltz, E. A. 2002. "Dark- 
SUSY: A numerical package for supersymmetric dark matter calculations." Preprint astro- 
ph/0211238. See |http:// www.physto.se/^edsjo/darksusyl 



Gondolo, P. & Gelmini, G. B. 1991. "Cosmic Abundances Of Stable Particles: Improved 
Analysis," Nucl. Phys. B360, 145 

Gondolo, P. & Silk, J. 1999. "Dark matter annihilation at the galactic center," Phys. Rev. 
Lett. 83, 1719 

Griest, K. & Seckel, D. 1991. "Three Exceptions In The Calculation Of Relic Abundances," 
Phys. Rev. D43, 3191 

Habig, A. [Super-Kamiokande Collaboration] 2001. "An indirect search for WIMPs with 
Super-Kamiokande." Preprint hep-ex/0106024 



48 



Hagiwara, K., et al. [Particle Data Group Collaboration] 2002. "Review Of Particle Physics," 
Phys. Rev. D66, 010001 

Hall, L. J., Moroi, T., & Murayama, H. 1998. "Sneutrino cold dark matter with lepton- 
number violation," Phys. Lett. B424, 305 

Hooper, D., & Dingus, B. L. 2002. "Limits on supersymmetric dark matter from EGRET 
observations of the galactic center region." Preprint astro-ph/0210617 

Ibata, R., Lewis, G. P., Irwin, M., Totten, E., & Quinn, T. 2001. "Great Circle Tidal 
Streams: Evidence for a Nearly Spherical Massive Dark Halo around the Milky Way," 
Astrophys. J. 551, 294 

Ibata, R. A., Wyse, R. F. G., Gilmore, G., Irwin, M. J., & Suntzeff, N. B. 1997. "The 
Kinematics, Orbit, and Survival of the Sagittarius Dwarf Spheroidal Galaxy," Astron. J. 
113, 634 

Jungman, G., Kamionkowski, M., & Griest, K. 1996. "Supersymmetric dark matter," Phys. 
Rep. 267, 195 

Kane, G. L., Wang, L. T., & Wells, J. D. 2002. "Supersymmetry and the positron excess in 
cosmic rays," Phys. Rev. D65, 057701 

Kim, J. E. 1979. "Weak Interaction Singlet And Strong CP Invariance," Phys. Rev. Lett. 
43, 103 

Kuzmin, V., & Tkachev, I. 1998. "Ultra-high energy cosmic rays, superheavy long-living 
particles, and matter creation after inflation," JETP Lett. 68, 271 

Lauer, T. R., Statler, T. S., Ryden, B. S., & Weinberg, D. H. 1986 (approx). "A New 
and Definitive Met a- Cosmology Theory." Scanned images available on Steve McGaugh's 
website at www.astro.umd.edu/~ssm/mond/flowchart.html. 

Ma, C.-P. 1999. "Neutrinos and Dark Matter." Preprint astro-ph/9904001 

Majewski, S. R., Skrutskie, M. P., Weinberg, M. D., & Ostheimer, J. C. 2003. "A Two 
Micron All Sky Survey View of the Sagittarius Dwarf Galaxy. I. Morphology of the Sagit- 
tarius Core and Tidal Arms," Astrophys. J. 599, 1082 

Markevitch, M., et al. 2003. "Direct constraints on the dark matter self-interaction cross- 
section from the merging galaxy cluster 1E0657-56." Preprint astro-ph/0309303 

Mayer- Hasselwander, H. A. et al. 1998. "High-energy gamma-ray emission from the Galactic 
Center," A&A 335, 161 

Merritt, D., Milosavljevic, M., Verde, L., & Jimenez, R. 2002. "Dark Matter Spikes and 
Annihilation Radiation from the Galactic Center," Phys. Rev. Lett. 88, 191301 

Miralda-Escude, J. 2002. "A Test of the CoUisional Dark Matter Hypothesis from Cluster 
Lensing," Astrophys. J. 564, 60 



49 



Moore, B., Governato, F., Quinn, T., Stadel, J., & Lake, G. 1998. "Astrophys. J. Lett. ," 
Resolving the Structure of Cold Dark Matter Halos499, L5 

Navarro, J. F., Frenk, C. S., & White, S. D. M. 1996. "The Structure of Cold Dark Matter 
Halos," Astrophys. J. 462, 563 

Newberg, H. J., Yanny, B., Grebel, E. K., Hennessy, G., Ivezic, Z., Martinez-Delgado, D., 
Odenkirchen, M., Rix, H.-W., Brinkmann, J., Lamb, D. Q., Schneider, D. P., & York, D. 
2003. "Sagittarius Tidal Debris 90 Kiloparsecs from the Galactic Center," Astrophys. J. 
Lett. 596, 191 

Peccei, R. D., & Quinn, H. R. 1977. "CP Conservation In The Presence Of Instantons," 
Phys. Rev. Lett. 38, 1440 

Persic, M., Salucci, P., & Stel, F. 1996. "The universal rotation curve of spiral galaxies - I. 
The dark matter connection," MNRAS 281, 27; Erratum zbzd. 283, 1102 (1996) 

Porquet, D., Predehl, P., Aschenbach, B., Grosso, N., Goldwurm, A., Goldoni, P., Warwick, 
R. S., & Decourchelle, A. 2003. "XMM-Newton observation of the brightest X-ray flare 
detected so far from Sgr A*," A&A 407, L17 

Primack, J. R. 2001. " Whatever Happened to Hot Dark Matter?," SLAC Beam Line31N3, 
50 [Preprint Eitro^h/0112336| 

Sarkar, S. 2003. "Neutrinos from the Big Bang." Preprint hep-ph/0302175 

Shifman, M. A., Vainshtein, A. I., & Zakharov, V. I. 1980. "Can Confinement Ensure Natural 
CP Invariance Of Strong Interactions?," Nucl. Phys. B166, 493 

Silk, J. 2003. Lectures at this School. 

Snowden-Ifft, D. P., Martoff, C. J., & Burwell, J. M. 2000. "Low pressure negative ion drift 
chamber for dark matter search," Phys. Rev. D61, 101301 

Spergel, D. N. & Steinhardt, P. J. 2000. "Observational Evidence for Self-Interacting Cold 
Dark Matter," Phys. Rev. Lett. 84, 3760 

Spergel, D. N., et al. 2003. "First- Year Wilkinson Microwave Anisotropy Probe (WMAP) 
Observations: Determination of Cosmological Parameters," Astrophys. J. Suppl. 148, 175 

Stiff, D., Widrow, L.M., & Frieman, J. 2001. "Signatures of hierarchical clustering in dark 
matter detection experiments," Phys. Rev. D64, 083516 

Tovey, D. R., Gaitskell, R. J., Gondolo, P., Ramachers, Y., & Roszkowski, L. 2000. "A new 
model-independent method for extracting spin-dependent cross section limits from dark 
matter searches," Phys. Lett. B488, 17 

Ullio, P., Bergstrom, L., Edsjo, J., & Lacey, C. 2002. "Cosmological dark matter annihila- 
tions into 7 rays: A closer look," Phys. Rev. D66, 123502 



50 



Verde, L., et al. 2002. "The 2dF Galaxy Redshift Survey: the bias of galaxies and the density 
of the Universe," MNRAS 335, 432 

Wald, R. M., 1994. "Quantum Field Theory In Curved Space-Time And Black Hole Ther- 
modynamics" (Chicago, IL: University of Chicago Press) 

Yoshida, N., Springel, V., White, S. D. M., & Tormen, G. 2000. "Weakly Self-interacting 
Dark Matter and the Structure of Dark Halos," Astrophys. J. Lett. 544, L87 

Zhitnitsky, A. R. 1980. "On possible suppression of the axion hadron interactions," Sov. J. 
Nucl. Phys. 31, 260 



51 



